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EXECUTIVE  SUMMARY 


A  series  of  open-air  dispersion  experiments  known  collectively  as  the 
Evaporation,  Recirculation,  and  Dispersion  in  Light  Winds  trials  was  conducted 
at  U.S.  Army  Dugway  Proving  Ground  (DPG)  in  July-August  1995.  The  partici¬ 
pants  were  scientists  representing  the  British  Chemical/Biological  Defence 
Establishment  (CBDE);  the  Canadian  Defence  Research  Establishment  Suffield 
(DRES) ;  U.S.  Army  DPG  West  Desert  Test  Center  (WDTC) ;  and  a  consortium  from 
the  University  of  Utah,  University  of  Maryland,  and  Michigan  State  University. 
These  researchers  examined  several  aspects  of  surface  effects  on  the  disper¬ 
sion^  of  windborne  material  within  the  lower  portion  of  the  atmospheric 
boundary  layer.  The  objectives  of  the  three  main  experimental  subtests  were 
to:  (1)  obtain  evaporation  data  from  neat  and  thickened  methyl  salycilate 

(MES)  and  diethyl  malonate  (DEM)  droplets  deposited  on  sand,  concrete,  and 
aluminum  surfaces;  (2)  docvunent  the  effects  of  atmospheric  stability  and 
obstacle  geometry  on  the  accumulation  and  dissipation  of  material  entrained 
into  the  recirculation  region  within  the  wake  of  a  2-m  cube;  and  (3)  examine 
iicsr-surface  turbulence  and  its  effects  on  dispersing  plvimes  during  low  wind 
speed  conditions .  Additional  methodology  studies  were  performed  to  determine 
the  efficacy  of  using  fast-response  hot-wire  anemometry  to  measure  the 
turbulent  wind  field,  including  turbulence  production  and  dissipation,  around 
obstacles  and  in  the  viscous  sublayer  immediately  adjacent  to  the  earth's 
surface.  Significant  results  were  obtained  from  all  of  the  subtests  and 
methodology  studies . 

The  Evaporation,  Recirculation,  and  Dispersion  in  Light  Winds  trials 
series  was  a  continuation  of  a  tripartite  (U.S.,  U.K. ,  and  Canada)  concentra¬ 
tion  fluctuation  test  series  designed  to  improve  our  understanding  of  atmo¬ 
spheric  dispersion  processes.  As  noted  above,  this  field  test  consisted  of 
three  subtests  and  two  methodology  studies.  The  Evaporation,  Recirculation, 
and  Plume  Profile  and  Turbulence  Effects  Subtests  were  conducted  at  the  DPG 
Tower  Grid  test  site.  The  methodology  studies  were  conducted  at  Tower  Grid 
and  a  remote  test  site  adjacent  to  the  Goodyear  Road  causeway,  which  crosses 
the  salt  flats  of  the  Great  Salt  Lake  Desert.  Figure  1  shows  the  locations  of 
test  sites  on  a  DPG  map.  The  data  from  the  Evaporation  and  Recirculation 
Subtests  were  analyzed  by  CBDE  contractors  who  have  documented  the  results  in 
reports  to  CBDE  and  a  doctoral  thesis.  The  methodology  study  results  have 
been  documented  in  University  of  Utah  technical  reports,  open  literature 
publications,  and  a  doctoral  thesis.  The  principal  focus  of  this  report  is  on 
the  Plume  Profile  and  Turbulence  Effects  Subtest,  which  studied  near- surface 
plume  dispersion  under  non-steady  flow  in  light  winds. 

Table  1  in  Section  2  summarizes  the  trials  conducted  during  the  Evapora¬ 
tion  Subtest.  Of  the  34  successful  trials  of  30-  to  120-min  duration,  12  used 
MES  as  the  simulant  and  22  used  DEM.  The  subtest  objective,  which  was  to 
document  the  effects  of  droplet  size,  thickener  content,  simulant  type, 
surface  t3rpe,  and  meteorological  conditions  on  droplet  evaporation  rate,  was 
met.  Data  from  these  trials  in  combination  with  data  from  earlier  work  have 
been  used  to  validate  a  comprehensive  evaporation  model,  satisfying  a  knowl¬ 
edge  shortfall  identified  by  Technical  Panel  9  of  The  Technical  Cooperation 
Program  (TTCP)  Subgroup  E  on  Chemical/Biological  (CB)  Defense.  Detailed 
results  are  provided  by  Roberts  (1996). 
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Figure  1.  Map  of  Dugway  Proving  Ground  showing  the  locations  of  the  Tower 


Table  2  in  Section  2  summarizes  the  89  trials  of  the  Recirculation  Region 
Subtest.  Forty-nine  of  these  trials,  which  were  tjrpically  2  to  5  min  in 
duration,  were  performed  with  a  cube  orientation  of  0®  (cube  face  normal  to 
the  wind)  and  44  were  performed  with  a  cube  orientation  of  45®.  Data  acquired 
during  the  trials  included  infrared  (IR)  imagery  of  tracer  gas  (propylene) 
dispersion  around  the  cube  and  hot-wire  anemometer  measurements  of  winds  and 
turbulence  within  the  cube's  wake.  The  subtest  objective,  which  was  to 
document  obstacle  geometry  and  atmospheric  stability  effects  on  the  retention 
of  windborne  material  entrained  into  the  accumulation  region  in  the  wake  of  a 
bluff  object,  was  met.  Flow  around  bluff  bodies  during  stable  atmospheric 
conditions  produces  large  wakes  that  retain  high  concentrations  of  windborne 
contaminants.  However,  when  these  results  are  normalized  by  wake  width  and 
wind  speed,  atmospheric  stability  appears  to  have  a  minimal  effect  on  contami¬ 
nant  residence  time.  Recirculation  trials  results  are  being  used  to  develop 
and  validate  urban  dispersion  model  algorithms.  Detailed  analyses  of  obstacle 
wake  entrainment  results  are  given  by  Mavoridis  (1996,  1997). 

The  university  consortium  conducted  several  methodology  studies  using  hot¬ 
wire  anemometry.  One  study,  which  supported  the  Recirculation  Region  Subtest, 
involved  turbulence  measurements  within  the  wake  of  the  2-m  cube.  This  study 
was  only  partly  successful  because  dust  accumulated  on  the  hot-wire  sensor 
elements.  The  second  methodology  study,  which  investigated  stability  and 
Re3molds  number  effects  on  flow  within  the  near -surface  wall  region,  was 
performed  on  the  salt  flats  immediately  north  of  the  Causeway  near  Photo  Pad 
11.  Steady  wind  conditions  and  the  absence  of  micrometer- scale  windborne  dust 
particles  made  high  quality  near -surface  hot-wire  anemometry  possible  at  this 
location.  These  studies  included  wind  velocity  profile  measurements  around  a 
small  (0.25-m)  cube  and  near-surface  turbulence,  vorticity,  and  stress 
measurements.  The  resulting  data  document  surface  (near-wall)  effects  on 
turbulence  and  momentum  profiles,  fluxes,  and  turbulent  kinetic  energy 
generation  and  dissipation  in  high  Reynolds  number  flows  and  provide  a 
database  for  development  of  near-surface  high  Reynolds  number  flow  models. 

The  Plume  Profile  and  Turbulence  Effects  Sub test  was  designed  to  provide 
information  on  the  effects  of  near -surface  turbulence  on  the  vertical  compo¬ 
nent  of  plume  dispersion  during  non-steady  flow  in  light  winds.  During  this 
subtest,  propylene  was  continuously  released  for  periods  of  15  to  35  min  at 
heights  ranging  from  0.05  to  2.5  m  above  ground  level  (AGL)  .  A  vertical  array 
of  gas  concentration  detectors  was  positioned  12,5  to  100  m  downwind  of  the 
source  to  provide  detailed  vertical  plume  concentration  profile  measurements. 
Concentration  detector  availability  limited  the  vertical  array  to  a  single 
tower  fully  instrumented  from  0.25  to  8.0  m.  Sonic  anemometer/thermometers 
(sonics)  also  mounted  on  this  tower  provided  wind  and  turbulence  measurements. 


Wind  meander  wafted  the  disseminated  propylene  plume  intermittently  across 
the  detectors  during  the  Plume  Profile  and  Turbulence  Effects  Subtest.  The 
quantity  of  nonzero  plume  concentration  measurements  varied  inversely  with  the 
amount  of  meander  experienced  during  each  trial.  Of  the  24  trials,  15 
provided  sufficient  propylene  concentration  data  to  form  vertical  concen¬ 
tration  profiles,  six  of  which  supported  detailed  statistical  analyses.  The 
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plume  vertical  profile  and  accompanying  micrometeorological  data  provide 
detailed  information  on  near-surface  dispersion  during  light  winds. 

The  Plume  Profile  and  Turbulence  Effects  Subtest  led  to  the  development  of 
two  new  indices,  the  Diabatic  Ratio  (DR)  and  the  Total  Diabatic  Influence 
(TDI),  DR  represents  the  ratio  of  convectively- driven  to  shear-driven 
turbulence,  and  TDI  is  a  measure  of  turbulence  strength.  Both  indices  are 
needed  to  describe  the  turbulent  state  of  the  atmosphere:  DR  defines  whether 
turbulence  is  primarily  convective  or  shear -driven,  while  TDI  defines  its 
magnitude.  These  indices  were  developed  as  alternatives  to  existing  similari¬ 
ty-based  indices  to  better  represent  the  state  of  the  atmosphere  within  the 
first  few  meters  of  the  surface  where  the  preponderance  of  troop  exposure  to 
windborne  hazards  occurs.  When  applied  to  realistic  non-steady  atmospheric 
conditions  such  as  light  winds  punctuated  by  turbulence  bursts  and  transition 
through  sunrise,  these  indices  appear  to  characterize  the  state  of  the  surface 
layer  better  than  the  widely  used  similarity  theory  indices.  DR  and  TDI  also 
can  be  easily  calculated  in  near  real-time  using  the  output  of  a  two -axis 
sonic  anemometer/thermometer .  If  these  new  indices  are  proven  useful  in 
further  testing,  they  could  find  wide  application  in  the  modeling  of  windborne 
hazards . 
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SECTION  1 .  INTRODUCTION 


The  Surface  Effects  on  Evaporation,  Recirculation,  and  Dispersion  in  Light 
Winds  project  was  designed  to  address  issues  related  to  troop  vulnerability  to 
high  concentrations  of  hazardous  windborne  materials .  The  principal  focus  of 
this  report  is  on  the  effects  of  surface  characteristics  (e.g.,  wakes,  poros¬ 
ity,  surface  heating,  and  turbulence)  on  the  concentrations  of  material 
dispersing  near  the  surface  during  light  wind  conditions.  Near- surface 
effects  are  important  because  the  majority  of  troop  exposures  to  hazardous 
materials  occur  within  a  few  meters  of  the  surface.  Field  measurements  are 
needed  to  define  these  effects  because  near -surface  velocity  gradients  are  too 
strong  and  turbulence  scales  are  too  small  for  present  direct  numerical  or 
large  eddy  simulations.  Surface  effects  are  a  significant  operational  concern 
to  the  military  because  high  concentrations  of  hazardous  materials  can  occur 
near  the  surface,  particularly  during  light  wind  situations  when  the  turbu¬ 
lence  driving  the  dispersion  process  is  weak.  A  major  goal  of  this  work  is  to 
improve  our  understanding  of  near-surface  dispersion,  which  should  lead  to 
improved  hazard  scenario  modeling  and  better  design  and  deployment  of 
windborne  hazard  monitoring  equipment  (e.g.,  chemical  and  biological  agent 
detectors) . 

The  Surface  Effects  on  Evaporation,  Recirculation,  and  Dispersion  in  Light 
Winds  experiment  series  was  a  tripartite  effort  with  participation  by  the 
British  Chemical/Biological  Defence  Establishment-Porton  (CBDE) ,  the  Canadian 
Defence  Research  Establishment  Suf field  (DRES),  and  U.S.  Army  Dugway  Proving 
Ground  (DPG) .  Each  participating  nation  contributed  unique  capabilities  to 
address  multiple  aspects  of  near-surface  dispersion.  The  principal  CBDE  focus 
was  on  evaporation  and  bluff  body  recirculation  region  studies,  while  DRES 
perfoinned  detailed  near-surface  concentration  profile  measurements.  DPG  was 
responsible  for  test  program  management  and  meteorological  and  infrared  (IR) 
imagery  support.  A  DPG  contractor,  the  University  of  Utah  Department  of 
Mechanical  Engineering,  provided  the  hot-wire  anemometer  measurements  needed 
to  characterize  near-surface  velocity  and  turbulence  profiles. 

The  Evaporation,  Recirculation,  and  Dispersion  trials  were  conducted  at 
two  locations,  one  near  Tower  Grid  (40°  06'N,  112°  58'W)  and  one  near  Photo 
Pad  11  on  the  salt  flats  immediately  north  of  Goodyear  Road  (40°  08 'N,  113° 
26'W).  Tower  Grid  is  typical  of  high  desert  terrain,  with  long  periods  of 
light  winds  persisting  through  the  evening  and  early  morning  hours.  The 
terrain  is  characterized  by  a  gently  sloping  clay  surface  with  a  light  cover 
of  short,  brushy  vegetation.  Drainage  flows  make  Tower  Grid  a  suitable 
location  for  dispersion  and  evaporation  tests  as  long  as  measurements  are  not 
required  immediately  above  the  dusty  clay  surface.  In  contrast  to  the 
aerodynamically  rough  Tower  Grid  location,  the  salt  crust  surface  on  the  Great 
Salt  Lake  Desert  forms  an  aerodynamically  smooth  surface  virtually  devoid  of 
vegetation.  Velocity  and  turbulence  profile  measurements  very  close  to  the 
surface  are  more  successful  on  the  salt  flats  because  the  salt  crust  is  free 
of  the  fine  clay  particles  that  interfere  with  the  calibration  of  the  fine 
heated  wires  used  as  sensing  elements  for  hot-wire  anemometers. 

The  desert  floor  forms  an  impenetrable  "no-slip"  boundary  underlaying  the 
atmospheric  boundary  layer  (ABL) ,  the  portion  of  the  troposphere  sufficiently 
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close  to  the  earth  that  it  responds  to  surface  forcings.  This  surface  is  a 
momentum  sink,  a  heat  source  during  the  day,  and  a  heat  sink  at  night. 

Because  the  constantly  varying  fluxes  of  heat  and  momentum  keep  ABL  flow  in  a 
continuous  state  of  evolution,  a  "steady  state"  condition  that  can  be  ade¬ 
quately  treated  using  standard  statistical  methodologies  is  the  exception 
rather  than  the  rule.  The  structure  of  the  ABL  can  be  divided  into  an  outer 
region,  an  inertial  sublayer,  and  a  near-wall  region.  Each  region  or  sublayer 
has  unique  characteristics  that  can  be  described  in  terms  of  relevant  scaling 
variables.  The  focus  of  this  report  is  on  the  inertial  sublayer,  which 
encompasses  the  first  few  tens  of  meters  above  the  surface,  and  the  near-wall 
region,  which  lies  within  a  few  tens  of  centimeters  of  the  surface. 

The  state  of  the  ABL  is  usually  described  in  terms  of  similarity  theory 
relationships,  with  Monin-Obukhov  scaling  applied  to  the  surface  layer,  while 
mixed  layer  or  convective  scaling  is  used  when  convection  is  present.  These 
semi -empirical  relationships  are  based  on  an  apparent  behavioral  dependence  of 
boundary  layer  characteristics  on  scaling  variables  that  include  the  vertical 
fluxes  of  heat  and  momentum.  Because  a  flux  describes  the  passage  of  a 
quantity  (heat,  momentum)  through  a  surface,  it  should  be  representative  of  an 
area.  In  contrast,  surface -based  measurements  used  to  calculate  fluxes  are 
usually  made  at  a  point  in  space,  which  renders  vertical  flux  calculations 
highly  susceptible  to  non- stationary  flows  over  inhomogeneous  surfaces.  This 
problem  is  particularly  prevalent  close  to  the  surface  where  strong  gradients 
and  constantly  changing  fluxes  of  heat  and  momentum  occur.  Thus,  while  many 
authors  (for  example,  Deardorff,  1985;  Venkatram,  1988)  have  successfully 
applied  similarity  scaling  at  "smoke  stack"  heights  of  10  to  100  m  and  through 
the  convection- dominated  boundary  layer,  no  comparably  successful  methodology 
has  emerged  for  the  light  winds  and  variable  conditions  that  frequently  occur 
at  "nose  height"  (within  2  m)  of  the  surface. 

Persistent  weak  winds  within  the  inertial  sublayer  present  a  dispersion 
modeling  challenge  because  some  of  the  greatest  hazards  for  windborne  releases 
of  toxic  materials  occur  during  light  winds  when  mixing  is  weak  or  intermit¬ 
tent  and  flow  patterns  are  variable.  High  concentrations  of  toxic  materials 
can  also  persist  for  long  periods  or  accumulate  in  the  wakes  of  obstacles 
under  these  conditions  in  complex  terrain.  Variable  light  winds  are  most 
common  during  near -surface  nocturnal  drainage  flows,  but  can  also  extend 
through  sunrise  transition  into  periods  of  daytime  convection.  This  report 
examines  the  limitations  of  present  theory  and  presents  some  practical 
alternatives  for  describing  the  state  of  the  inertial  sublayer  during  variable 
light  wind  episodes . 

The  near-wall  region  of  the  ABL  has  received  little  attention  by  the 
meteorological  community  in  part  because  of  the  difficulty  of  obtaining 
measurements  close  to  the  earth's  surface.  However,  ABL  turbulence  is 
intimately  coupled  with  surface  heat  and  momentum  fluxes.  Also,  fine-scale 
processes  such  as  evaporation,  deposition,  resuspension,  and  the  destruction 
of  concentration  gradients  occur  as  a  consequence  of  near-wall  turbulence 
generation  and  dissipation.  This  report  describes  flow  and  turbulence 
measurements  within  the  near -wall  region  and  the  use  of  the  resulting  data  to 
verify  near -surface  physical  models. 
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SECTION  2.  THE  EVAPORATION  AND  RECIRCULATION  MEASUREMENT  PROGRAM 


2.1  THE  EVAPORATION  SUBTEST 

The  Evaporation  Subtest  was  undertaken  to  provide  a  validation  data  base 
for  a  mass  transfer  model  for  predicting  the  evaporation  of  chemical  agent 
droplets  from  various  surfaces  (Roberts,  1994).  This  model  considers  the 
effects  on  droplet  evaporation  of  surface  roughness,  molecular  diffusion 
within  the  wall  region  adjacent  to  the  surface,  and  atmospheric  stability. 

The  development  and  validation  of  this  model  helps  to  fill  a  modeling  short¬ 
fall  identified  by  an  expert  panel  of  Technical  Panel  9  of  The  Technical 
Cooperation  Program  (TTCP)  Subgroup  E  on  Chemical/Biological  (CB)  Defense. 

The  need  for  a  new  evaporation  model  was  driven  by  the  unsatisfactory 
performance  of  existing  models.  Roberts  (1994)  reports  that  available  droplet 
evaporation  results  can  be  reproduced  with  existing  models  only  through  the 
use  of  unrealistic  spread  factors,  wind  velocity  adjustments,  and/or  absorp¬ 
tion  or  diffusion  conditions  because  the  model  physics  is  inadequate. 
Evaporation  model  development  has  been  stymied  by  incomplete  characterization 
of  the  absorption  and  evaporation  processes ,  including  the  lack  of  an  adequate 
model  of  near- surface  turbulence. 

The  specific  objectives  of  the  Evaporation  Subtest  were  to:  (a)  quantify 
spread  factors  and  evaporation  rates  for  droplets  deposited  on  various 
surfaces,  and  (b)  obtain  a  better  understanding  of  the  relative  importance  of 
various  physical  factors  (droplet  size,  surface  porosity,  simulant  type,  and 
thickener  content)  on  these  variables.  The  surfaces  chosen  for  the  trials 
were  sand,  concrete,  and  aluminum  plate,  which  respectively  represented 
porous,  semi-porous,  and  non-porous  substrates.  Diethyl  Malonate  (DEM,  a 
simulant  for  the  chemical  agent  GD)  and  Methyl  Salicylate  (MES,  a  simulant  for 
the  chemical  agent  HD)  were  the  selected  simulants.  Polymethylmethacrylate 
powder  (co-pol3mier  K-125)  served  as  the  thickener.  The  mean  droplet  size 
varied  from  1.0  to  3.6  mm,  depending  on  thickener  content. 

The  evaporation  trials  were  conducted  at  the  DPG  Tower  Grid  test  site 
between  17  July  and  11  August  1995.  The  site,  30  m  south  of  the  Tower  Grid 
Command  Post  (CP)  Building,  consisted  of  a  concrete  annular  ring  "sidewalk”  of 
6-m  radius  and  1-m  width  located  on  flat,  open  terrain  characterized  by  a 
surface  roughness  length  z^  of  1.6  cm.  The  interior  of  the  annulus  was  filled 
with  sand,  and  aluminum  plates  ringed  its  outside  edge.  The  purpose  of  the 
annular  ring  design  was  to  ensure  uniform  exposure  of  detectors  stationed  at 
the  annulus  centroid  to  the  evaporated  vapors  from  droplets  deposited  in  a 
uniform  pattern  around  the  ring.  Ultraviolet  ion  collectors  were  placed  at 
the  annulus  center  at  heights  of  5,  15,  30,  50,  and  80  cm  above  ground  level 
(AGL)  to  measure  these  evaporated  vapors.  Ancillary  measurements  obtained 
near  the  evaporation  site  included  surface  and  air  temperature  profiles.  A 
sonic  anemometer/thermometer  (sonic)  mounted  at  2  m  AGL  measured  wind  and 
turbulence.  A  schematic  of  the  annular  ring  and  experimental  set-up  is  shown 
in  Figure  2. 

DEM  or  MES  admixed  with  0,  3,  or  5  percent  by  weight  of  co-pol3rmer  K125 
thickener  were  deposited  in  the  trials  series  summarized  in  Table  1.  The 
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Table  1.  Evaporation  Site  Trial  Sxiiranary.* 


Trial 

No. 

Date 

Day/Mo/Year 

Duration 

(min) 

Droplet 

Diameter 

(mm) 

Liquid 

Type 

Liquid 
Loading** 
(mg  m  *) 

Thickener 

(%) 

Surface 

Type 

uOI 

21/7/95 

60 

3.6 

DEM 

3358 

0 

concrete 

u02 

21/7/95 

60 

3.6 

DEM 

7750 

0 

concrete 

u03 

22/7/95 

90 

3.6 

DEM 

7750 

0 

concrete 

u04 

22/7/95 

90 

0.2-2.0 

DEM 

9376 

3 

concrete 

u05 

24/7/95 

60 

02.-2.0 

DEM 

9055 

3 

aluminium 

u06 

24/7/95 

45 

0.2-2.0 

DEM 

6781 

5 

aluminium 

u07 

24/7/95 

30 

3.6 

DEM 

4421 

0 

aluminium 

u08 

24/7/95 

30 

3.6 

DEM 

4421 

0 

aluminium 

u09 

24/7/95 

30 

0.2-2.0 

DEM 

6781 

5 

aluminium 

ulO 

24/7/95 

30 

0.2-2.0 

DEM 

7127 

3 

aluminium 

Ul1 

25/7/95 

60 

3.6 

DEM 

7233 

0 

concrete 

u12 

25/7/95 

90 

0.2-2.0 

DEM 

9035 

3 

concrete 

u14 

25/7/95 

60 

3.6 

DEM 

6006 

0 

concrete 

u15 

26/7/95 

90 

0.2-2.0 

DEM 

8261 

3 

sand 

u16 

26/7/95 

60 

0.2-2.0 

DEM 

4721 

5 

sand 

u18 

28/7/95 

60 

0.2-2.0 

DEM 

6119 

3 

sand 

u19 

28/7/95 

60 

0.2-2.0 

DEM 

5993 

5 

sand 

u20 

28/7/95 

60 

3.6 

DEM 

7327 

0 

sand 

u22 

29/7/95 

90 

0.2-2.0 

DEM 

5071 

5 

concrete 

u23 

29/7/95 

60 

0.2-2.0 

DEM 

5178 

3 

concrete 

u24 

29/7/95 

60 

0.2-2.0 

DEM 

5993 

5 

sand 

u25 

29/7/95 

60 

0.2-2.0 

DEM 

6884 

3 

sand 

u26 

01/8/95 

90 

0.2-2.0 

MES 

4734 

3 

concrete 

u27 

01/8/95 

120 

0.2-2.0 

MES 

4636 

5 

concrete 

u28 

01/8/95 

90 

3.6 

MES 

6928 

0 

concrete 

u30 

02/8/95 

60 

0.2-2.0 

MES 

5680 

3 

sand 

u31 

02/8/95 

75 

0.2-2.0 

MES 

6823 

5 

sand 

u32 

03/8/95 

75 

0.2-2.0 

MES 

5679 

3 

aluminium 

u35 

03/8/95 

60 

0.2-2.0 

MES 

8036 

5 

aluminium 

u36 

04/8/95 

90 

3.6 

MES 

5298 

0 

aluminium 

u37 

04/8/95 

90 

3.6 

MES 

9211 

0 

sand 

u38 

04/8/95 

60 

0.2-2.0 

MES 

10061 

3 

sand 

u39 

04/8/95 

45 

0.2-2.0 

MES 

7883 

5 

sand 

u40 

04/8/95 

45 

0.2-2.0 

MES 

6509 

3 

concrete 

•  From  Roberts  (1996). 

Liquid  ioading  is  the  neat  chemical  loading  per  unit  length  of  annulus. 
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droplet  deposition  was  fixed  at  10  g/m^  with  droplet  diameters  ranging  from 
0.2  to  2  mm  (mean  of  1.0  mm)  for  thickened  mixtures  and  3.6  mm  for  neat  mix¬ 
tures.  The  micro-pipette  that  deposited  the  neat  simulant  droplets  provided 
greater  droplet  size  control  than  the  rotating-disk  droplet  sprayer  used  with 
the  thickened  mixture.  A  uniform  evaporating  (curved)  line  source  was  formed 
by  depositing  a  continuous  swath  of  droplets  at  a  fixed  rate  along  the  arc, 
which  provided  uninterrupted  detector  exposure  to  a  uniform  evaporating  source 
even  with  changes  in  wind  direction.  A  video  recorder  provided  time -lapsed 
images  of  deposited  droplet  samples  as  they  soaked  into  and/or  evaporated  from 
their  underlying  surfaces. 


Figure  2.  Schematic  of  the  annular  ring  and  experimental  set-up  used  in  the 
Evaporation  Trials  (from  Roberts,  1996). 
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The  concentrations  of  the  evaporating  simulants  were  measured  by  CBDE's 
Ultraviolet  Ion  Collector  (UVIC)  photo- ionization  detectors  (PIDs).  As  part 
of  data  analysis,  the  UVIC  voltage  readings  were  converted  to  simulant 
concentrations  in  parts  per  million  (ppm)  and  the  horizontal  vapor  flux  was 
calculated  for  each  trial.  The  vapor  flux  was  computed  by  integrating  the 
product  of  concentration  and  wind  speed  over  the  plume  depth  using  a  power- law 
wind  profile  fitted  to  the  measured  wind  speeds  and  a  reflected  Gaussian 
distribution  fitted  to  the  measured  simulant  concentrations.  This  vapor  flux 
represented  the  mass  release  rate  per  unit  line  source  length  of  the  evaporat¬ 
ing  droplets.  Assuming  that  individual  droplets  contribute  equally  to  the 
total  amount  of  evaporating  material ,  the  calculated  flux  and  known  number  of 
individual  droplets  per  unit  length  of  line  source  define  the  droplet  evapora¬ 
tion  for  each  trial.  These  calculated  evaporation  rates  were  also  correlated 
with  time-lapse  droplet  size  video  imagery. 

The  quality  of  the  evaporation  data  was  evaluated  by  comparing  the  mass  of 
simulant  initially  deposited  on  the  aluminum  surfaces  to  the  vapor  mass 
measured  throughout  each  trial.  The  average  recovery  for  drops  deposited  on 
the  aluminum  surface  was  96  percent,  with  a  standard  deviation  of  15  percent. 
Roberts  (1996)  believes  that  similar  recovery  rates  were  achieved  for  drops 
deposited  on  the  sand  and  concrete  surfaces.  The  major  source  of  uncertainty 
occurred  toward  the  end  of  each  test  when  the  vapor  concentrations  were  very 
low. 

Roberts  (1996)  compares  the  Evaporation  Sub test  evaporation  data  to  the 
predictions  of  the  Roberts  (1994)  model  and  finds  that  the  model  reasonably 
replicates  the  early  stages  of  evaporation,  correctly  distinguishing  the 
effects  of  thickener  and  surface  porosity  on  evaporation  rates.  There  is  a 
close  (within  ±10  percent)  correspondence  between  predictions  and  measurements 
in  the  early  stages  of  evaporation,  with  perhaps  a  slight  early  stage 
underprediction  of  evaporation  on  concrete  surfaces.  However,  the  model  tends 
to  overpredict  the  evaporation  rate  in  the  later  stages  of  evaporation 
regardless  of  the  surface  type  or  thickener  content.  Roberts  (1996)  believes 
that  the  later  stages  of  evaporation  are  governed  by  transport  processes 
within  the  surface  substrate  (such  as  capillarity,  diffusion,  and  gravitation¬ 
al  flow)  ,  and  that  the  physics  of  these  processes  is  not  well  specified.  More 
detailed  models  and  measurement  techniques  are  needed  to  resolve  this  issue. 

The  Evaporation  Subtest  results  also  provided  an  opportunity  to  quantify 
droplet  spread  in  terms  of  a  spread  factor,  the  ratio  of  the  wet  spot  diameter 
on  the  surface  to  the  droplet  diameter  prior  to  impact.  Droplet  images  prior 
to  and  after  impact  were  obtained  using  a  video  camera.  Unfortunately,  poor 
contrast  over  the  aluminum  surface  precluded  differentiating  wet  spots  on  that 
surface.  However,  Roberts  (1996)  was  able  to  achieve  the  required  video  image 
pairs  over  sand  and  concrete.  He  reports  a  spread  factor  of  2.5  on  sand 
regardless  of  droplet  size  or  viscosity.  In  contrast,  viscosity  and  droplet 
size  appear  to  affect  droplet  size  on  concrete.  A  rough  estimate  of  the 
spread  factor  on  concrete  is  0.6d  +  2.5,  where  d  is  the  droplet  diameter  in 
millimeters . 
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2.2  THE  RECIRCULATION  SUBTEST 


Flow  around  a  large  obstacle  forms  a  horizontally-oriented  "horseshoe" 
vortex  that  begins  upwind  of  the  obstacle  front  and  wraps  around  the  sides, 
trailing  into  the  wake.  This  horseshoe  vortex  mixes  windborne  contaminants 
rapidly  to  the  ground,  distributing  them  around  the  obstacle  and  into  a 
recirculating  "accumulation  region."  Contaminant  concentration  within  the 
accumulation  region  varies  as  a  function  of  wake  dimensions,  wind  speed, 
turbulence  intensity,  and  the  magnitude  and  duration  of  the  contaminant 
concentration  impinging  upon  the  obstacle.  Gaussian  dispersion  models  do  not 
presently  include  methodologies  to  account  for  the  accumulation,  retention, 
and  dissipation  of  contaminants  trapped  within  the  wakes  of  large  obstacles. 

To  address  this  deficiency,  Mavoridis  (1996,  1997)  performed  a  series  of  wind 
tunnel  experiments  and  field  trials,  including  the  Recirculation  Subtest.  The 
principal  purpose  of  the  Recirculation  Subtest  was  to  investigate  source 
geometry  and  atmospheric  stability  effects  on  the  accumulation  and  dissipation 
of  a  windborne  tracer  (propylene)  in  the  accumulation  region  within  the  wake 
of  a  2-m  pl3rwood  cube. 

The  Recirculation  Subtest  was  conducted  at  the  DPG  Tower  Grid  test  site 
between  22  July  and  8  August  1995.  Trials  were  conducted  on  flat,  open 
terrain  characterized  by  low  brush  and  a  surface  roughness  length  of  1.6  cm. 
The  test  site  was  approximately  800  m  south  of  the  Evaporation  Subtest  site. 
This  separation  was  sufficient  to  preclude  the  possibility  that  tracer 
materials  released  at  one  site  would  contaminate  measurements  made  at  the 
other. 

Each  Recirculation  Subtest  trial  began  with  the  positioning  of  the 
disseminator  4  or  5  m  upwind  of  the  cube  at  a  height  of  1  or  2  m  AGL.  As 
shown  in  Figure  3,  the  cube,  which  was  mounted  on  a  central  pivot,  could  be 
rotated  to  present  either  a  "face-on"  (0°  orientation)  or  "corner-on"  (45® 
orientation)  front  to  the  mean  wind  direction.  UVICs  were  stationed  in  the 
wake  of  the  cube  at  a  height  of  1  m  to  obtain  time -re solved  measurements  of 
propylene  concentration.  Table  2  summarizes  the  Recirculation  Subtest  trials. 
The  last  trial  set,  designated  UTH,  in  Table  2,  is  distinguished  by  the 
location  of  the  source  0.5  m  behind  the  cube  as  shown  in  Figure  4.  This 
configuration  produced  reference  maximum  concentrations  within  the  accumula¬ 
tion  region. 

A  recirculation  trial  could  only  be  accomplished  while  the  tracer  source 
remained  aligned  upwind  of  the  cube  front.  Thus,  trial  conduct  required 
careful  monitoring  of  the  meandering  winds.  Once  the  wind  direction  came  into 
alignment,  the  disseminator  was  turned  on  to  begin  a  trial.  The  UVICs  were 
then  switched  on  to  sample  the  propylene  concentrations  within  the  accumula¬ 
tion  region.  As  this  propylene  concentration  reached  its  upper  plateau,  the 
disseminator  was  abruptly  shut  off.  The  UVICs  then  measured  the  propylene 
concentration  decay  within  the  accumulation  region.  The  multiple  trials 
conducted  within  experimental  sessions  included  several  source  distances , 
source  heights,  and  cube  orientations. 

The  recirculation  region  propylene  concentration  measurements  were  supple¬ 
mented  by  sonic  anemometer  readings  obtained  from  a  CEDE -furnished  tripod- 
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mounted  sonic  positioned  upwind  of  the  cube  at  a  height  of  2-m  AGL.  The  sonic 
provided  mean  wind  and  turbulence  statistics  and  the  fluxes  of  heat  and 
momentum  required  for  determinatin  of  the  Obukhov  length  L.  Recirculation 
Subtest  micrometeorological  summaries  are  provided  by  Mavoridis  (1996,  1997). 

A  total  of  89  trials  of  2-  to  5 -min  duration  were  completed  during  the 
Recirculation  Subtest.  Concentration  data  from  UVICs  stationed  in  the  cube 
wake  were  used  to  estimate  the  accumulation  region  residence  time  t^,  which  is 
defined  by  Mavoridis  (1996)  as  the  time  required  for  the  accumulation  region 
concentration  to  decay  to  1/e  (37  percent)  of  its  orginal  value.  The  trial 
results  Indicate  that  accumulation  region  evacuation  proceeds  exponentially, 
as  predicted  by  Hunt  and  Castro  (1994) .  Figure  5  illustrates  an  averaged 
recirculation  region  time  series  with  this  characteristic  exponential  decay 
rate. 

The  Recirculation  Subtest  residence  time  t^  results,  when  normalized  by  the 
cube  dimension  h  and  mean  flow  u  around  the  cube,  were  initially  found  to  be 
inconsistent  with  results  presented  by  Fackrell  (1984),  particularly  when  the 
cube  was  in  a  side -on  orientation  with  respect  to  the  wind  direction. 

However,  the  hot-wire  anemometer  measurements  described  by  Metzger  and 
Klewicki  (1996,  see  also  Section  4  of  this  report)  showed  that  the  side-on 
cube  orientation  produced  a  larger  wake  width  h'  than  the  face -on  orientation. 
When  Mavoridis  (1996)  defined  a  non-dimensional  residence  time  t^  as 


tj.  =  t^u/h'  . 


(2-1) 


he  was  able  to  achieve  results  consistent  with  Fackrell  (1984).  Thus,  h' 
appears  to  be  a  more  relevant  dimension  than  h  for  wake  accumulation  modeling. 

Flow  around  the  cube  during  stable  atmospheric  conditions  produced  larger 
wakes,  higher  accumulation  region  concentrations,  and  longer  t^  than  corre¬ 
sponding  measurements  taken  when  the  atmosphere  was  unstable.  However,  the 
magnitude  of  t^  varied  little  with  stability,  indicating  that  u  and  h' 
accounted  for  most  of  the  variability  observed  in  t^.  Mavoridis  (1996) 
concludes  that  stability  effects  are  implicitly  accounted  for  in  suitably 
normalized  accumulation  region  residence  times. 
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Figure  3.  (a)  Plan  view  of  Recirculation  Subtest  source,  UVICs,  and  2-in  cube, 

with  cube  orientations  of  0®  (face-on,  solid  line)  and  45°  (side- 
-on,  dashed  line),  (b)  Side  view  of  Recirculation  Subtest  source, 
UVICs,  and  2-m  cube  with  cube  orientation  of  0°. 
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Figure  4.  Plan  view  of  source  and  detector  configuration  for  Experimental 

Session  UTM  conducted  on  7  August  1996  with  the  source  positioned 
at  a  1.0  m  height  in  the  center  rear  of  the  2-m  cube  and  an  array 
of  UVICs  stationed  within  the  accumulation  region  in  the  wake  of 
the  cube,  (b)  Side  view  of  configuration  for  Experimental  Session 
UTM. 
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Time  (s) 


Figure  5.  Averaged  concentration  time  series  from  UVICs  stationed  in  a  cross- 
wind  array  1.5  m  downwind  of  the  2-m  cube  during  Trial  UTE02, 
illustrating  exponential  decay  following  gas  switch  off  at  1063  s 
into  the  trial  sequence  (from  Mavoridis,  1996). 


Table  2.  Recirculation  Subtest  Trial  Summary. 
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SECTION  3,  ATMOSPHERIC  BOUNDARY  LAYER  DYNAMICS 


3.1  BOUNDARY  LAYER  THEORY  AND  PRACTICES 

The  atmospheric  boundary  layer  (ABL)  is  the  portion  of  the  atmosphere 
characterized  by  the  high  Reynolds  number  flows  that  are  produced  by  the 
diurnal  heat  flux  cycle  and  other  forcings  originating  at  the  earth's  surface. 
The  ABL  is  generally  more  complex  than  "engineering"  boundary  layers  (found  in 
pipe,  channel,  and  wind  tunnel  flows)  because  it  contains  continuously  varying 
heat  fluxes,  temperature  and  velocity  gradients,  directional  shear  and 
divergence  fields,  and  complex  surface  roughness  features.  Perhaps  the  ABL's 
major  distinguishing  feature  is  that  it  is  in  a  constant  state  of  change, 
continuously  adjusting  to  varying  flxixes  of  heat  and  momentum.  Also,  the 
spatial  and  temporal  scales  parameterized  by  the  Re)molds  niomber  relevant  to 
ABL  flow  are  often  several  orders  of  magnitude  greater  than  comparable  scales 
found  in  other  boundary  layer  flows.  ABL  theory  is  further  complicated  by  the 
lack  of  a  unified  treatment;  the  outer  region  and  inertial  sublayer  are 
typically  studied  by  meteorologists  interested  in  atmospheric  transport  and 
diffusion,  while  the  wall  region  remains  the  domain  of  mechanical  engineers 
who  are  principally  interested  in  momentum  and  turbulence  exchange  in  highly 
sheared  flows  and  viscous  effects  on  near- surface  interactions.  This  section 
summarizes  salient  features  of  existing  boundary  layer  theories  and  practices, 
followed  by  a  discussion  of  their  limitations  under  non-steady  (variable) 
light  wind  conditions. 

Figure  6  is  a  schematic  illustration  of  the  ABL  evolution  during  a  diurnal 
cycle  over  high  desert  terrain.  The  principal  ABL  regions  are  depicted  as 
they  evolve  in  space  and  time.  Wavy  lines  at  region  boundaries  represent 
fluctuations  in  boundary  heights  due  to  perturbations  such  as  gravity  waves 
and  convection.  The  principal  ABL  scaling  variables  (defined  later  in  this 
section)  are  shown  in  brackets .  The  dashed  and  dotted  lines  denote  boundaries 
between  stable,  unstable,  and  free  convection  zones.  The  depth  and  duration 
of  the  various  regions  and  layers  depicted  in  Figure  6  vary  with  factors  such 
as  heat  flux,  surface  roughness,  and  the  pressure -driven  velocity  field. 

The  state  of  the  ABL  is  too  complex  and  poorly  understood  to  be  derived 
from  first  principles.  Consequently,  current  ABL  theory  is  based  on  empirical 
relationships  and  Reynolds  approximations  to  the  Navier- Stokes  equations  that 
allow  the  use  of  certain  similarity  theory  "scaling  variables"  to  describe 
salient  characteristics  of  the  ABL.  Pertinent  scaling  variables  can  be  formed 
into  dimensionless  parameters  that  merge  scattered  data  into  "self-similar" 
forms  (universal  curves,  functions,  or  constants)  that  depict  ABL  attributes. 
The  principal  scaling  variables  and  dimensionless  parameters  for  the  outer 
region,  inertial  sublayer,  and  wall  region  are  described  below  using  the  usual 
micrometeorological  Cartesian  convention  in  which  downwind  and  crosswind 
direction  distances  are  represented  respectively  by  x  and  y,  with  z  depicting 
height  above  ground  level.  Wind  velocity  components  in  the  alongwind, 
crosswind,  and  vertical  directions  are  represented  by  u,  v,  and  w,  respective- 

ly. 
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Figure  6 .  Schematic  illustration  of  ABL  evolution  during  a  diurnal  cycle 
over  high  desert  terrain.  The  principal  ABL  scaling  variables 
are  shown  in  brackets.  Dashes  and  asterisks  denote  boundaries 
between  stable,  unstable,  and  free  convection  zones. 
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Boundary  layer  theory  is  developed  under  the  following  assumptions:  (1) 
flow  can  be  represented  as  a  2 -dimensional  velocity  field  with  a  uniform  mean 
^®^tical  gradient;  (2)  velocity,  temperature,  and  pressure  fluctuations  are 
statistically  homogeneous;  (3)  all  statistical  quantities  for  meteorological 
variables  are  stationary  with  respect  to  time;  and  (4)  viscous  dissipation 
outside  the  wall  region  is  small.  The  stationarity  assvunption,  wherein  all  of 
the  a/at  terms  of  the  following  equations  are  zero,  is  particularly  problemat¬ 
ic  in  a  boundary  layer  characterized  by  constant  change.  Given  these  assump¬ 
tions  and  neglecting  terms  higher  than  second  order,  simplified  Reynolds 
approximations  for  turbulent  shear  stress  {-WiP)  and  vertical  heat  flux 
(w'T')  budgets  are  given  respectively  by  (Garratt,  1992) 


3i?i7/3t  =  -^du/dz  -  3uV^/3z  +  {g/T)  u'T'  -  p"^  (  u'dP^/dz  +  vPdP'/dx)  (3-1) 

and 

dvFr’/dt  =  -W^dd/dz  -  dw^T^/dz  +  (g/T)!^  -  (T^dP'/dz)  (3-2) 


where  6  is  the  potential  temperature,  g  is  gravitational  acceleration,  T  is 
absolute  temperature,  [1]  is  a  shear  (du/dz)  or  stratification  (dB/dz) 
production  term,  [2]  is  the  vertical  turbulent  flux  divergence  of  uV  or  iTr^, 
[3]  represents  buoyancy  production  (conversion  of  turbulent  kinetic  energy  to 
turbulent  potential  energy) ,  and  [4]  is  a  pressure  interaction  term  which 
tends  to  dissipate  shear  stress  and  heat  flux.  Overbars  and  primes  denote 
time  averages  and  departures  therefrom,  respectively.  Applying  similar 
treatment  to  the  three  velocity  variance  components  produces  the  turbulent 

kinetic  energy  =  (u^  +  +  w^)/2.  The  resulting  conservation  equation  is 


dg^/dt  =  -u'vPdu/dz  +  {g/D  w^T'  -  d{w'cf  +  vPp'/p)  /bz  -  e  (3-3) 


where  [1]  represents  shear  production,  [2]  depicts  buoyant  production  or 
dissipation,  [3]  consists  of  vertical  turbulent  energy  transport  and  pressure- 
velocity  correlation  terms,  and  [4]  is  the  turbulent  kinetic  energy  dissipa¬ 
tion  term  e .  Most  scaling  variables  and  parameters  used  to  describe  the  state 
of  the  ABL  relate  to  the  production,  flux  divergence,  and  buoyancy  terms  of 
Equations  (3-1)  through  (3-3). 

3.2  THE  OUTER  REGION 

The  outer  region  is  the  upper  portion  of  the  ABL.  Its  top  extends  from  a 
few  tens  of  meters  AGL  during  quiescent  nocturnal  conditions  to  several 
kilometers  AGL  during  strong  convection.  Horizontal  winds  in  the  outer  region 
respond  principally  to  the  interaction  of  pressure  gradient  and  Coriolis 
forces,  with  the  surface  influence  mainly  expressed  as  convection  initiated  by 
surface  heating.  Exchange  between  the  outer  region  and  other  regions  within 
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the  ABL  continues  so  long  as  convection  is  sufficiently  vigorous  to  distribute 
heat,  momentum,  and  other  dispersing  scalar  quantities  between  the  surface  and 
the  top  of  the  mixed  layer.  The  top  of  :  he  mixed  layer  Zi  is  usually  defined 
by  a  capping  inversion.  The  presence  oi  .onvective  mixing  is  the  principal 
feature  that  distinguishes  the  outer  region  from  the  free  atmosphere  above  it. 
Wind  velocity  profiles  within  the  outer  region  are  independent  of  surface - 
based  wind  shear. 

The  outer  region  expands  during  the  day  as  the  heat  flux  H  from  the 
earth's  surface  produces  convective  motions  that  raise  the  capping  inversion. 
The  sign  and  magnitude  of  H  is  related  to  the  vertical  velocity- temperature 
covariance  Vr'  by 


H  =  pCpVT>, 


(3-4) 


where  p  is  atmospheric  density  and  Cp  is  the  specific  heat  of  air  at  constant 
pressure.  A  positive  w'T'  indicates  an  upward  flux  of  heat  from  the  surface. 
When  a  sufficient  positive  heat  flux  is  present  to  maintain  convection,  the 
outer  region  becomes  a  convection- dominated  "mixed  layer."  Mixed  layer 
similarity  relationships,  which  have  fairly  uniform  vertical  velocity  variance 
profiles  and  minimal  surface  drag  effects,  predominate  between  about  O.lZi  and 

0. 8zi. 

Heat  flux  and  z^  are  the  most  important  scaling  variables  for  the  ABL  outer 
region.  They  can  be  combined  to  form  the  convective  velocity  scale 


gz^H 


PC^T 


1/3 


(3-5) 


which  characterizes  the  mixed  layer  state  of  the  outer  region.  The  convective 
velocity  scale  is  tjrpically  zero  at  sunrise,  increases  rapidly  to  a  mid¬ 
afternoon  maximum  between  1.5  and  3.0  m/s,  and  decreases  toward  zero  in  the 
evening  (Stull,  1988,  Figure  4.1).  The  significance  of  w.  is  its  relevance  to 
the  scaling  of  outer  region  velocity  variables.  For  example,  the  vertical 
velocity  variance  is  a  measure  of  the  vigor  with  which  vertical  mixing 
occurs  within  the  ABL.  Although  direct  measurements  of  in  the  outer 
region  are  rarely  available,  the  convective  contributions  to  ^  normalized  by 
w,^  fall  onto  a  universal  curve  of  height  z  normalized  by  Zi  according  to  the 
expression  (Stull,  1988) 


Thus,  if  Zi  and  v/T'  are  known.  Equation  (3-6)  provides  an  analytical  expres¬ 
sion  for  the  vertical  variation  of  in  the  convective  ABL. 
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When  the  ratio  of  the  mean  wind  speed  U  to  the  convective  velocity  scale  w, 
falls  within  the  range  2<U/w.<10,  dispersion  within  the  upper  regions  of  the 
ABL  can  be  parameterized  using  mixed- layer  scaling.  The  time  scale  for  large 
convective  eddies  is  z^/w. ,  and  the  "travel  time"  for  dispersing  material  is 
rendered  dimensionless  through  multiplication  by  w./z^.  For  long  range 
dispersion,  downwind  travel  distance  x  from  a  dissemination  position  scales 
with  the  product  of  w,/Zi  and  x/U,  forming  a  dimensionless  downwind  distance 
(Deardorff,  1985) 

y  nr 

X  =  ^  .  (3-7) 

UZi 


When  plotted  against  X,  plume  dispersion  within  the  convective  ABL  exhibits 
consistent  universal  behaviors.  For  example,  laboratory  and  field  studies 
(Kaimal  et  al.,  1986)  indicate  that  the  centroid  of  dispersing  material 
released  at  the  surface  into  a  convective  ABL  tends  to  rise  within  X  <  1.0  of 
the  release  point,  reaching  a  height  of  0.8z/Zi  between  1.0<X<2.0.  It 
then  gradually  descends  toward  0.5  z/Zi  at  X  >  3.0,  becoming  well  mixed  at 
greater  downwind  distances . 

In  the  absence  of  convection,  the  boundary  layer  depth  is  sustained  only 
by  mechanical  turbulence  generated  by  flow  over  roughness  elements.  As 
surface  heating  and  convection  subside  in  the  evening,  the  outer  region 
collapses  into  layers  isolated  from  each  other  by  temperature  inversions.  The 
nocturnal  mixing  height  is  within  the  lowest  of  these  inversions.  Shear 
layers  develop  when  turbulent  mixing  becomes  weak.  Above  the  nocturnal  mixing 
height,  the  atmosphere  slumps  into  a  "residual  layer"  (Stull,  1988)  character¬ 
ized  by  multiple  shear  zones  and  inversions  which  sandwich  layers  of  air  that 
are  decoupled  from  the  surface,  each  other,  and  the  free  atmosphere  above 
them.  Local  similarity  (Wyngaard,  1973)  is  dominant  within  each  residual 
layer.  Diffusing  material  released  within  one  of  these  layers  can  travel  a 
considerable  distance  without  mixing  into  the  layer  above  or  below  it. 

3.3  THE  INERTIAL  SUBLAYER 

Between  the  convection-dominated  mixed  layer  and  the  viscosity- dominated 
wall  region  lies  the  inertial  sublayer.  This  layer,  also  known  as  the  surface 
layer  or  the  constant  flux  layer,  is  tjTpically  a  few  meters  to  a  few  tens  of 
meters  in  depth.  The  inertial  sublayer  depth  is  the  closest  analog  to  the 
thickness  of  a  laboratory  boundary  layer  because  it  is  the  depth  through  which 
surface -generated  shear  forces  affect  flow  within  the  ABL. 

The  inertial  sublayer  is  so  named  because  it  supports  an  inertial  sub¬ 
range,  a  band  of  turbulent  eddy  scales  ranging  from  the  integral  scale  \j  to 
the  Taylor  microscale  Xx.  The  inertial  subrange  includes  eddy  scales  small 
enough  to  decay  isotropically,  distributing  energy  evenly  in  all  dimensions, 
but  large  enough  to  be  unaffected  by  viscosity.  Consequently,  inertial 
subrange  eddy  motions  are  the  same  order  of  magnitude  in  all  directions . 

The  integral  scale  defines  eddies  on  the  large  end  of  the  inertial 
subrange.  These  eddies,  with  a  characteristic  length  of  tens  to  hundreds  of 
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meters,  extract  kinetic  energy  from  the  mean  flow.  They  contain  the  greatest 
turbulent  kinetic  energy  and  are  believed  to  be  primarily  responsible  for  the 
turbulent  mixing  of  diffusing  material.  The  breakup  of  integral  scale-sized 
eddies  initiates  a  transfer  of  turbulent  energy  down  the  range  of  eddy  scales 
xmtil  they  are  dissipated  in  molecular  motions  as  heat. 

Integral  scale  length  depends  on  factors  such  as  distance  from  the  surface 
and  atmospheric  stability.  A  first  approximation  of  X,  is  given  by  Prandtl's 
hypothesized  characteristic  eddy  length  scale  of  kz  (Schlicting,  1979),  where 
k  is  the  von  Karman  constant  (0.4  ±0.02)  and  z  is  height  AGL.  The  largest 
integral  scales  are  associated  with  convective  eddies  in  the  upper  part  of  the 
inertial  sublayer. 

Near  the  lower  end  of  the  inertial  subrange  lies  the  Taylor  microscale, 
where  eddies  are  small  enough  (on  the  order  of  centimeters)  to  be  affected  by 
viscosity.  The  Taylor  microscale  is  defined  in  terms  of  q^,  £,  and  kinematic 
viscosity  v  as 

A,j.  =  (vg2/e) .  (3-8) 


A  continuous  cascade  of  energy  occurs  between  Xj  and  X^  as  large  eddies  break 
up  and  distribute  their  energy  isotropically  over  progressively  smaller 
scales.  Scales  smaller  than  X,.  form  a  gradual  transition  from  the  inertial 
subrange  to  the  dissipation  subrange,  where  viscous  effects  predominate  and 
the  destruction  of  concentration  gradients  by  molecular  diffusion  occurs . 
Similarity  arguments  work  well  in  the  inertial  sublayer  when  a  stationary 
inertial  subrange  extends  between  the  integral  scale  and  Taylor  microscale. 

The  inertial  sublayer  is  most  often  characterized  using  Monin-Obukhov 
similarity  theory,  which  assumes  that  the  average  values  of  buoyancy  and  the 
fluxes  of  heat  and  momentum  are  constant.  The  principal  Monin-Obukhov 
variables  are  the  heat  flxix  H,  roughness  length  Zj,,  friction  velocity  u., 
height  z,  and  ratio  of  gravitational  acceleration  to  absolute  temperature  g/T. 
The  effects  of  surface  roughness  are  parameterized  by  Zg.  Roughness  length 
accounts  for  the  effects  of  surface  texture  and  roughness  element  spacing  on 
the  inertial  sublayer  wind  profile.  The  friction  velocity,  which  is  the 
square  root  of  the  vertical  momentum  flux 

u*  =  V  -u'l/  ,  (3-9) 


is  the  most  important  Monin-Obukhov  velocity  scale. 

The  Obukhov  length  L,  formed  from  the  ratio  of  the  buoyancy  and  shear - 
driven  contributions  to  turbulent  kinetic  energy  (terms  [1]  and  [2]  in 
Equation  (3-3)),  is  the  principal  Monin-Obukhov  length  scale.  The  ratio  of  z 
to  L  forms  the  nondimens ional  stability  parameter 

The  magnitude  of  z/L  is  determined  primarily  by  the  ratio  of  heat  and  momentum 
fluxes,  while  heat  flux  determines  the  sign  of  z/L. 
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(3-10) 


_z  _  -  kg  zH 
^  pCpTul 

Monin-Obukhov  similarity  theory  describes  inertial  sublayer  mean  flow 
properties  as  functions  only  of  z/L.  For  -z/L  <1,  inertial  sublayer  velocity 
fluctuations  are  on  the  order  of  tu  (typically  0.1  to  1.0  m/s)  and  the  wind 
velocity  increases  linearly  with  the  logarithm  of  height  normalized  by  Ze. 

The  normalized  wind  speed  profile  is  expressed  as 


where  V'm  is  a  z/L- dependent  diabatic  influence  function.  Equation  (3-11) 
remains  valid  through  the  inertial  sublayer  and  into  the  wall  region  until  the 
flow  is  severely  distorted  by  near-surface  viscous  effects. 

Another  useful  Monin-Obukhov  relationship  is  the  ratio  of  the  vertical 
velocity  standard  deviation  to  u, .  The  average  a„  normalized  by  u.  is 
nearly  constant,  exhibiting  a  weak  z/L  dependence  that  becomes  important  only 
in  the  presence  of  strong  convection.  The  normalized  mean  profile  is  given 
by 


“1.3  +  f{z/L) 


(3-12) 


It  must  be  emphasized  that  Equation  (3-12)  represents  an  average  condition  and 
that  any  single  realization  of  the  ff«/u,  ratio  can  exhibit  a  deviation  of  ±50 
percent  or  more  due  primarily  to  the  inherent  temporal  variability  of 
An  analysis  performed  by  Wyngaard  (1973)  concludes  that  requires  samp- 

ling  periods  10  to  100  times  longer  than  those  for  the  vertical  velocity 
variance  to  achieve  comparable  levels  of  statistical  stability.  Consequently, 
Equation  (3-12)  is  an  attractive  alternative  to  the  eddy  correlation  method 
for  the  calculation  of  u.. 

The  dimensionless  ratio  of  height  (either  z  or  Zi)  to  L  is  widely  used  as  a 
measure  of  the  state  of  the  ABL  and  its  capacity  for  dispersion.  When  heat 
flux  is  weak,  L  increases  towards  infinity  and  the  atmosphere  approaches  an 
adiabatic  stratification  in  which  there  is  no  net  gain  or  loss  in  buoyancy  as 
air  parcels  rise  or  sink.  Motions  in  an  adiabatic  ABL  are  shear-driven,  being 
neither  enhanced  by  convection  nor  suppressed  by  density  stratification. 
Turbulent  energy  in  an  adiabatic  atmosphere  is  extracted  from  large  scale 
motions  during  flow  over  roughness  elements.  Diabatic  conditions  exist  when 
the  heat  flux  is  nonzero ;  the  magnitude  of  L  decreases  in  proportion  to  the 
magnitude  of  H.  Upward  heat  flux  produces  a  negative  L  associated  with 
convection- enhanced  dispersion.  Downward- directed  heat  flux  produces  a 
positive  L  with  damped  (negatively-buoyant)  turbulent  motions  and  suppressed 
vertical  dispersion.  Wyngaard  and  Cote'  (1971)  suggest  that  the  removal  of 
energy  by  damping  in  the  stable  surface  layer  is  locally  compensated  by 
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turbulent  energy  generated  near  the  surface  and  imported  into  the  inertial 
sublayer  by  the  turbulence  transport  term  [3]  in  Equation  (3-3).  Conversely, 
the  turbulence  transport  term  transports  energy  toward  the  surface  during 
convection. 

A  positive  (upward)  heat  flux  in  the  inertial  sublayer  quickly  grows 
through  the  outer  layer,  simultaneously  diminishing  -L  and  increasing  z^. 
Deardorff  (1972)  defines  a  -Zi/L  of  4.5  as  the  point  of  transition  between  a 
shear -driven  ABL  and  a  convection- driven  one  in  which  an  outer  region  mixed 
layer  overlies  a  convectively  unstable  inertial  sublayer.  When  the  heat  flxix 
is  sufficiently  strong  to  diminish  z/L  to  -1,  u'  and  w'  become  decorrelated 
and  u,  vanishes  in  free  convection.  Dispersing  material  released  into  free 
convection  typically  detaches  quickly  from  the  surface  and  diffuses  aloft  in 
vigorous  mixing.  A  free  convection  velocity  scale  Wf,  analogous  to  the  mixed 
layer  convective  velocity  scale  w. ,  is  the  pertinent  velocity  scale  in  an 
inertial  sublayer  dominated  by  free  convection.  Holtslag  and  Nieuwstadt 
(1986)  define  Wf  as 


= 


gzH 

pq,r 


11/3 


(3-13) 


Turbulent  motions  originating  in  the  free  convection  layer  blend  seamlessly 
into  the  mixed  layer  above  it. 


Parameters  independent  of  Monin-Obukhov  similarity  also  describe  the  state 
of  the  inertial  sublayer.  These  parameters  are  typically  formed  from  ratios 
of  fluxes,  velocity  and  temperature  variances  (aj,  ,  and  gradients  of 

velocity  (3u/3z)  and  potential  temperature  (dO/dz).  The  most  important  of 
these  non- similarity  parameters  include  the  Flux  Richardson  number,  eddy 
viscosity,  turbulent  thermal  diffusivity,  and  velocity  and  temperature 
correlation  terms.  The  Flux  Richardson  Number  Rp,  the  ratio  of  turbulent 
energy  production  by  density  variations  to  turbulent  energy  production  by  mean 
shear,  is  given  by  (Tavoularis  and  Corrsin,  1981) 


;  =  l£\ 

^  \  ^/  -u^v/du/dz 


(3-14) 


The  turbulent  or  eddy  viscosity,  which  represents  the  strength  of  turbulent 
diffusion  by  inertial  subrange -scale  eddies,  is  given  by 
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-u^w' 

du/dz 


(3-15) 


The  most  significant  terms  from  the  thermal  diffusivity  tensor  are  Dji  and  D22, 
which  represent  the  important  horizontal  and  vertical  components  of  turbulent 
thermal  diffusivity.  These  are  given  by 
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(3-16a) 
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(3-16b) 


The  velocity  and  temperature  correlation  terms,  which  describe  the  efficiency 
with  which  the  vertical  exchanges  of  heat  and  momentum  take  place,  are 
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v/t' 
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(3-17a) 
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(3-17b) 
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(3-17c) 


Neither  Monin-Obukhov  similarity  theory  nor  the  alternative  parameteriza- 
tions  work  well  outside  a  limited  range  of  conditions.  In  particular,  the 
strongly  stable  ABL  lacks  a  fully  developed  inertial  subrange  link  between 
large-scale  energetic  motions  and  the  Taylor  microscale  Thus,  the  cascade 

of  energy  down  the  eddy  scale  range  is  interrupted  in  the  strongly  stable  ABL. 
A  critical  flux  Richardson  number  defines  the  conditions  under  which  negative 
buoyancy  is  sufficient  to  break  the  inertial  subrange  link.  When  Rp  >  0.2 
there  is  insufficient  turbulence  to  maintain  a  continuous  transfer  of  heat  and 
momentum,  and  the  ABL  separates  into  strata  isolated  from  each  other  by  shear 
zones  and  inversions.  As  this  occurs,  z  ceases  to  be  a  relevant  scaling 
variable  and  a  state  of  "z-less  stratification"  (Wyngaard,  1973)  arises.  This 
condition,  with  vertical  mixing  confined  to  sporadic  turbulence  bursts, 
prevails  near  the  surface  during  nights  with  light  and  variable  winds.  In 
quiescent  periods  between  turbulence  bursts,  diffusing  material  released  into 
a  strongly  stable  layer  can  slump  to  the  surface  if  it  is  sufficiently  dense, 
or  persist  at  high  concentrations  confined  within  the  layer  into  which  it  is 
released. 

While  similarity  theory  descriptions  of  the  strongly  stable  ABL  have  long 
been  recognized  as  problematic,  it  is  increasingly  evident  that  Monin-Obukhov 
similarity  theory  assumptions  are  also  violated  during  periods  with  the  non¬ 
steady  conditions  associated  with  light  winds  regardless  of  stability.  Light 
winds  produce  great  uncertainty  in  flux  computations .  For  example ,  Biltoft 
(1993)  encountered  near-zero  or  counter -gradient  momenttim  fluxes  near  the 
surface  during  the  sunrise  transition  and  into  the  initial  development  of 
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convection.  When  uV  is  zero  or  positive,  u.  cannot  be  computed.  The 
problems  in  using  eddy  correlation  measurements  of  uV  to  estimate  u,  were 
first  recognized  by  Busch  and  Panofsky  (1968).  They  consider  eddy  correlation 
determinations  of  TPS?  to  be  unreliable  when  u.  is  less  than  0.32  m/s,  a  value 
that  is  rarely  reached  with  light  winds.  A  basic  problem  with  both  u't/  and 
is  that  they  are  the  algebraic  sums  of  positive  and  negative  contribu¬ 
tions  which,  in  many  circumstances,  substantially  cancel.  Thus,  these  fluxes 
can  be  small  compared  to  their  uncertainties  or  to  the  magnitudes  of  natural 
variations  arising  from  the  lack  of  homogeneity  and  stationarity .  Because  the 
magnitude  and  sign  of  L  are  largely  determined  by  the  ratio  of  u/  to  iPF, 
uncertainties  in  the  computation  or  estimation  of  those  variables  translate 
directly  into  uncertainties  in  L. 

The  familiar  ABL  parameterizations  not  associated  with  Monin-Obukhov 
similarity  theory  are  also  of  limited  use  during  periods  with  non-steady 
conditions.  For  example,  velocity  gradients  in  a  real  ABL  can  diminish  to 
zero  and  change  sign,  rendering  the  ABL  parameterizations  that  contain  them 
mathematically  intractible.  Wind  azimuth  and  elevation  angle  standard 
deviations  or  variances  and  wind  velocity  component  variances  frequently  are 
used  as  direct  measures  of  ABL  turbulence.  However,  because  these  variances 
are  subject  to  trends  and  meander,  it  is  possible  to  obtain  a  variance  of 
equivalent  magnitude  from  measurements  made  during  periods  of  low  wind  speed 
with  high  meander  and  measurements  made  during  periods  of  high  wind  speed  and 
high  turbulence  even  though  dispersion  under  these  two  regimes  will  be  very 
different.  Without  resorting  to  spectrxan  analysis,  it  is  not  possible  to 
determine  how  much  each  eddy  scale  within  the  flow  contributes  to  the  total 
variance .  Likewise ,  a  very  large  temperature  variance  can  be  measured  when  a 
steep  trend  or  abrupt  temperature  change  occurs,  even  when  there  is  little 
eddy  thermal  exchange.  Filtering  can  limit  trend  and  meander  contributions  to 
the  variance,  but  the  results  vary  as  a  function  of  the  technique  chosen  and 
filter  length.  Of  all  the  variances,  is  the  most  reliable  turbulence 
indicator  because  it  is  relatively  free  of  trends  and  large-scale  meander. 

Also,  the  surface-layer  similarity  relationship  between  and  the  rate  of 
vertical  dispersion  from  surface  or  near-surface  sources  is  well  established 
(see,  for  example,  Briggs,  1988;  Venkatram,  1988).  However,  vertical  velocity 
measurements  require  specialized  instruments  that  are  not  widely  available, 
and  vertical  motions  cannot  be  unambiguously  defined  in  uneven  or  sloping 
terrain. 

The  shortcomings  of  the  existing  inertial  sublayer  scaling  variables  and 
their  derived  parameterizations  suggest  a  need  for  more  robust  alternatives. 
Section  5  of  this  report  includes  a  detailed  examination  of  turblence  case 
studies  obtained  during  light  winds  and  proposes  alternative  methods  for 
parameterizing  near-surface  stability  and  turbulence  effects  on  diffusion. 

3. A  THE  WALL  REGION 

A  transition  to  the  wall  region  begins  at  the  bottom  of  the  inertial 
sublayer.  The  wall  region  is  composed  of  a  buffer  sublayer  and  a  viscous  or 
interfacial  layer  immediately  adjacent  to  the  surface.  Wind  profiles  remain 
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logarithmic  with  respect  to  height  into  the  buffer  sublayer,  but  the  vertical 
dimensions  of  turbulent  eddies  become  constrained  by  the  proximity  of  the 
surface  and  are  typically  smaller  than  eddy  horizontal  dimensions.  Also,  the 
effects  of  individual  roughness  elements  are  felt  within  the  buffer  sublayer, 
creating  very  inhomogeneous  turbulent  flow.  The  roughness  length  typically 
lies  within  the  buffer  sublayer  and  is  its  principal  length  scale.  Garratt 
(1994)  observes  that  the  transition  into  the  wall  region  t3rpically  begins  when 
z  <  lOOZo,  reaching  the  lower  limit  of  the  logarithmic  wind  profile  at 
approximately  lOZo. 

The  transition  from  the  inertial  subregion  to  the  buffer  sublayer  requires 
a  change  in  scaling  variables  as  viscosity  effects  become  dominant  while  z^ 
and  g/T  cease  to  be  relevant.  Because  little  information  exists  on  wall 
region  diabatic  effects,  the  relevance  of  the  heat  flxix  as  a  scaling  variable 
remains  largely  unknown.  The  composition  of  stress  also  changes  from  Reynolds 
stress  (u.)  to  surface  shear  stress  (u,) ,  the  product  of  the  velocity  gradient 
with  the  dynamic  viscosity  coefficient,  due  to  the  increased  contribution  from 
viscosity.  The  u./u,  ratio  approaches  unity  when  going  from  the  inertial 
sublayer  through  the  buffer  sublayer  as  the  contributions  to  u,  from  viscous 
stress  compensate  for  the  loss  of  Re3niolds  stress  near  the  surface.  Long  and 
Chen  (1981)  suggest  the  existence  of  a  u./u,  maximxam  at  a  distance  from  the 
surface  proportional  to  the  square  root  of  the  Reynolds  number,  placing  it 
near  the  boundary  between  the  buffer  sublayer  and  the  inertial  sublayer  in 
atmospheric  flows .  Although  the  existence  of  such  a  maximvim  is  not  yet 
rigorously  confirmed,  its  distance  from  the  surface  could  become  an  important 
length  scale.  Dimensional  quantities  measured  within  the  wall  region  are 
typically  rendered  nondimens ional  (inner-normalized)  using  u,  and/or  viscosity 

y .  The  vertical  distance  from  the  surface  can  be  represented  by  the  non- 
dimensional  height 

z*  =  zu^/v  .  (3-18) 


Using  inner  normalization,  the  buffer  sublayer  lies  within  5  <  z*  <  n,  where  n 
ranges  from  50  to  100.  The  viscous  sublayer  (z*  <  5)  forms  the  interface  with 
the  surface. 

The  viscosity  of  air  v  at  20  °C  is  0.000015  mVs,  which  is  four  or  five 
orders  of  magnitude  smaller  than  the  t3q)ical  eddy  viscosity  v,  in  the  inertial 
sublayer.  Therefore,  viscosity  effects  on  inertial  subrange -scale  eddies  are 
usually  ignored  outside  the  wall  region.  Garratt  (1992)  remarks  that  v^/v  is, 
like  the  Rejmolds  number,  a  ratio  of  inertial  to  viscous  forces  operating 
within  the  ABL.  However,  strong  gradients  within  the  wall  region  greatly 
diminish  v,  as  v  becomes  dominant.  Thus,  the  inertial  subrange  terminates  at 
the  Taylor  microscale  (on  the  order  of  0.1  m)  in  the  buffer  sublayer.  The 
dominant  eddies  in  the  viscous  sublayer  are  on  the  order  of  the  Kolmogorov 
scale  (0.001  m) . 

The  buffer  sublayer  is  close  enough  to  the  surface  that  flow  within  it 
experiences  strong  surface  roughness  effects ,  but  far  enough  from  the  surface 
that  it  does  not  experience  severe  viscous  retardation.  Therefore,  the 


23 


generation  of  the  small  scale  turbulence  responsible  for  the  rapid  breakdown 
of  concentration  gradients  reach's  its  maximum  in  the  buffer  sublayer. 

The  viscous  sublayer  lies  inm  sdiately  above  the  surface.  Logarithmic  mean 
wind  profiling  fails  in  the  viscous  sublayer  because  of  the  effects  of 
viscosity  and  the  presence  of  extreme  gradients  above  an  impermeable  surface. 
Motion  in  the  viscous  sublayer  only  occurs  when  the  momentum  is  sufficient  to 
overcome  the  resistance  threshold  of  retarding  viscous  forces.  Thus,  the  most 
prevalent  motions  within  the  viscous  sublayer  tend  to  be  intermittent  bursts 
aligned  along  the  axial  velocity. 

A  primary  consideration  for  near -wall  flow  is  whether  the  wall  is  aerody- 
namically  rough  or  smooth.  The  aerod3mamic  nature  of  the  surface  is  deter¬ 
mined  by  surface  roughness,  as  represented  by  texture  or  grain  size  (h,)  in  a 
ratio  with  the  surface  shear  stress  u^  and  viscosity  v.  If  Ujhj./v  >  75,  which 
is  the  usual  case  along  the  ground  surface,  the  flow  is  aerodynamically  rough 
with  turbulence  penetrating  to  the  surface  (Garratt,  1992).  In  this  case, 
both  small  scale  eddy  motions  and  molecular  diffusion  play  major  roles  in 
surface-to-air  heat  and  momentum  exchanges,  driving  processes  such  as  deposi¬ 
tion  and  evaporation.  The  Tower  Grid  site  is  aerod)mamically  rough,  but  the 
Salt  Flats  site  with  a  grain  size  of  1  mm  or  less  is  aerod3niamically  smooth 
(,Urh^/v  <  5)  under  most  conditions. 

Models  of  the  wall  region  (for  example,  Falco,  1991)  are  based  on  studies 
performed  using  fast-response  hot-wire  anemometers  in  shear-driven  wind  ttmnel 
flows  (see,  Willmarth  and  Lu,  1972).  The  primary  mechanism  for  energy 
exchange  is  the  turbulent  "burst"  in  which  "streaks"  of  low- speed  fluid  near 
the  surface  are  violently  "ejected"  into  the  lower  part  of  the  buffer  layer, 
to  be  displaced  by  downward  "sweeps"  of  higher  momentum  fluid.  The  downward 
moving  fluid  diverges  into  "pockets"  upon  contact  with  the  surface.  Pocket 
formation  forces  the  surrounding  low-momentum  air  to  converge  and  stretch  into 
thin,  elongated  "streaks"  with  vertical  acceleration,  re-initiating  the 
ejection  process.  Re3molds  stress  formation  occurs  along  the  upstream 
boundaries  of  the  ejected  low  speed  fluid  packets  as  they  enter  the  buffer 
sublayer,  mix  with  higher  speed  fluid,  and  turn  into  the  direction  of  the  mean 
flow.  Falco  (1991)  also  states  that  the  sweep/ejection  process  is  responsible 
for  a  majority  of  the  mechanical  turbulent  kinetic  energy  production  within 
the  wall  region. 

The  momentum  deficit  thickness  Reynolds  ntjunber  R,  describes  the  ratio  of 
inertial  forces  contributing  to  turbulent  motions  within  the  ABL  to  retarding 
viscous  forces.  This  Reynolds  number  is  given  by 

iJe  =  Ujijv  ,  (3-19) 


where  U.  is  the  free -stream  velocity  at  the  top  of  the  surface  layer  and  h^  is 
the  momentum  deficit  thickness.  Momentxan  deficit  thickness  is  an  "equivalent 
zero  velocity  distance"  used  as  an  indicator  of  the  amount  of  momentum  lost 
with  respect  to  the  free  stream  velocity  due  to  the  presence  of  the  boundary 
layer.  The  momentvim  deficit  thickness  is  roughly  equal  to  one -tenth  of  the 
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shear-driven  boundary  layer  thickness  (Klewicki  et  al.,  1995).  Using  this 
approximation,  R,  for  the  surface  layer  is  on  the  order  of  10‘. 
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SECTION  4.  WALL  REGION  STODTKS 


4.1  TEST  SITE  DESCRIPTION 

A  series  of  preliminary  near-surface  flow  studies  were  performed  by  a 
University  Consortium  (principally  members  of  the  University  of  Utah,  Univer¬ 
sity  of  Maryland,  and  Michigan  State  University  Mechanical  Engineering 
Departments)  at  DPG's  Surface  Layer  and  Environmental  Sciences  Test  (SLTEST) 
site,  which  is  1300  m  west  of  Photo  Pad  11  and  30  m  north  of  Goodyear  Road 
(see  Figure  1).  This  site  on  the  bottom  of  a  dry  lake  bed  was  selected  for 
surface  layer  studies  because  of  its  flat,  undisturbed,  vegetation- free 
surface  and  an  upwind  fetch  that  is  unobstructed  by  wake -generating  obstacles 
for  many  tens  of  kilometers.  Flow  over  the  dry  lake  bed  is  nearly  homoge¬ 
neous,  and  the  absence  of  wake -generating  obstacles  upwind  of  the  SLTEST  site 
minimizes  terrain- induced  meander  experienced  at  more  complex  sites.  Composed 
of  a  sand/clay/salt  mixture  capped  by  a  thin  salt  crust,  the  lake  bed  surface 
has  a  roughness  length  of  less  than  0.001  m  and  is  relatively  free  of 
micrometer -size  dust  particles  which  can  impinge  upon  and  alter  the  calibra¬ 
tion  of  micrometer- scale  hot-wire  probes.  A  rise  of  1  m  over  a  distance  of  10 
km  to  the  south  also  produces  a  negligible  (1/10000)  slope,  thereby  minimizing 
localized  gravity -driven  flows. 

In  addition  to  a  low  slope  angle ,  minimal  roughness ,  and  steady  winds ,  the 
SLTEST  site  is  characterized  by  large  Rj,  estimated  by  Metzger  and  Klewicki 
(1996)  to  be  on  the  order  of  several  million.  An  R*  of  this  magnitude  is 
representative  of  flows  within  the  ABL.  The  large  Rj  is  principally  a  conse¬ 
quence  of  the  large  h^  (on  the  order  of  tens  of  meters),  which  is  several 
orders  of  magnitude  greater  than  the  centimeter -scale  boundary  layers  generat¬ 
ed  within  typical  engineering  wind  tunnels.  Because  of  its  smooth  surface  and 
unobstructed  flow,  the  SLTEST  site  is  an  open  air  equivalent  of  a  large  wind 
tunnel  where  high  R,  flows  of  geophysical  interest  can  be  measured  without  the 
need  for  extremely  high  velocities  or  resorting  to  other  unrealistic  flow 
scaling  criteria. 

4.2  SLTEST  STUDIES 

The  University  Consortium  conducted  measurement  campaigns  at  the  SLTEST 
site  on  24-26  May  and  8  August  1994  and  on  21  July- 9  August  and  13  September 
1995.  Five  distinct  sets  of  experiments  were  conducted  during  these  cam¬ 
paigns:  (1)  viscous  sublayer  velocity  and  pressure  measurements,  (2)  axial 

and  vertical  velocity  profile  and  gradient  measurements  within  the  near-wall 
region,  (3)  velocity  and  pressure  measurements  around  a  cube,  (4)  flow 
visualization  of  coherent  structures  within  the  viscous  sublayer,  and  (5) 
studies  of  near-surface  turbulence  generation  and  dissipation. 

The  hot-wire  sensors  were  copper-plated  tungsten  wires  with  a  total 
length  of  3  mm  and  an  unplated  "active  region"  of  5-pm  diameter  and  1-mm 
length  maintained  at  a  constant  overheat  ratio  set  between  1.3  to  1.7.  Each 
probe  consisted  of  one  or  more  wires,  with  multi-wire  probes  designed  to 
sample  multiple  velocity  components  and/or  velocity  gradients  and  vorticity. 
When  velocity  profiles  were  needed,  the  hot-wire  probes  were  mounted  on  stands 
in  "rakes"  at  heights  ranging  from  0.5  mm  to  2.0  m  AGL.  The  hot-wire  systems 
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were  operated  at  sampling  rates  varying  between  500  and  10,000  Hz  for  trial 
durations  of  2  to  20  min.  Klewicki  et  al.  (1995)  and  Folz  (1997)  provide  a 
more  detailed  discussion  of  the  inrtrumentation  and  data  collection.  On-site 
hot-wire  calibrations  were  performed  immediately  before  and  after  each  trial 
using  a  pitot-static  tube  velocity  standard  and  differential  pressure  trans¬ 
ducer  housed  in  a  portable  planar  jet  facility.  The  planar  jet  facility  was 
built  by  the  University  of  Utah  Department  of  Mechanical  Engineering  specifi¬ 
cally  for  operational  use  in  the  open  atmosphere . 

The  flow  visualization  measurement  station,  which  is  shown  in  plan  view  in 
Figure  7,  consisted  of  a  polyethylene  slab  (1.3  x  2.0  x  0.025  m)  mounted  flush 
with  the  salt  flat  surface.  The  slab  contained  an  injection  slit  through 
which  theatrical  fog  can  be  disseminated  to  form  a  uniform  curtain  immediately 
above  the  slab  surface.  A  video  camera  recorded  the  behavior  of  this  thin 
smoke  curtain  while  pressure  measurements  were  taken  with  Larson-Davis  Model 
2250  0.25 -in  microphones  embedded  within  the  polyethylene  slab.  These  micro¬ 
phones  respond  to  pressure  fluctuations  at  frequencies  in  the  range  of  23  Hz 
to  100  kHz.  Thus,  the  flow  visualization  measurement  station  provides  the 
capability  to  characterize  activity  within  the  viscous  sublayer. 
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Figure  7.  Plan  view  schematic  of  the  SLTEST  site  experimental  setup,  includ¬ 
ing  the  polyethylene  slab,  smoke  generation  system,  and  stand  sup¬ 
porting  the  hot-wire  probes  and  camera  (from  Klewicki,  et  al., 
1995). 


4.3  WALL  REGION  TRIAL  RESULTS 


Direct  axial  velocity  gradient  and  shear  stress  measurements  were  made  by 
aligning  hot-wire  probe  arrays  into  the  mean  flow  direction  within  the  viscous 
sublayer  and  buffer  layer.  The  acquired  data  were  compiled  into  statistical 
summaries  that  included  trial  means,  root-mean- square  (rms)  deviations  about 
the  mean,  skewness,  and  kurtosis.  All  results  were  presented  as  a  function  of 
z*,  and  compared  with  previous  results  obtained  in  lower  R,  wind  tunnel  flows. 

The  inner-normalized  mean  axial  velocity  (u*  -  u/u,)  within  the  viscous 
sublayer  was  found  to  increase  linearly  with  height,  followed  by  a  transition 
to  a  logarithmic  profile  beginning  at  z*  =  10  (see  Figure  8) .  Also  shown  in 
Figure  8  is  a  solid  line  representing  Coles  "law  of  the  wall"  for  logarithmic 
wind  profiles  in  the  buffer  sublayer,  which  is  given  by  (Murlis  et  al.,  1982) 
as 

u*  =  -^  =  5.61  Inz*  +  f(i?e)  (4-1) 


where  f(R«)  is  weakly  dependent  on  the  Reynolds  ntjunber  (*.5.0  according  to 
Klewicki  and  Falco,  1990).  The  data  from  this  study  (Folz,  1997;  Metzger  and 
Klewicki,  1996)  depart  significantly  from  this  velocity  profile  model,  but  are 
in  reasonable  agreement  with  the  dotted  lines  in  Figure  8,  which  represent  the 
results  presented  by  Blackwelder  and  Haritonidis  (1983). 

The  rms  velocity  variations  found  within  the  viscous  sublayer  in  this 
study  are  virtually  equivalent  to  those  found  in  previous  measurements,  but 
the  buffer  layer  rms  values  (see  Figure  8b),  are  nearly  33  percent  larger  than 
those  measured  previously  in  lower  R^  flows.  Likewise,  the  skewness  and 
kurtosis  are  proportionally  larger,  particularly  within  the  viscous  sublayer, 
than  previously  measured  (see  Figures  8c  and  d)  .  Klewicki  et  al.  (1995) 
postulate  that  the  increased  buffer  layer  rms  velocity  (u')  may  indicate  an  Rj 
dependence  in  inner-normalized  axial  velocity  rms  (u'*  «=  u'/u,)  and/or  im¬ 
proved  probe  spatial  resolution  and  surface  shear  stress  definition.  They 
also  suggest  that  the  strong  positive  skewness  in  u'  is  a  consequence  of  the 
no-slip  condition  existing  at  the  surface. 

Klewicki  et  al.  (1995)  describe  how  the  ejection  slit  smoke  curtain  was 
used  to  produce  visual  images  of  "streaks"  and  "pockets"  forming  within  the 
viscous  sublayer  flow.  The  smoke  released  in  a  continuous,  uniform  thin 
curtain  at  the  surface  quickly  converged  into  long,  stretched  filaments 
(streaks)  of  low  speed  fluid  bordered  by  U-shaped  clear  pockets  of  clear  fluid 
recently  swept  into  the  viscous  sublayer  from  above.  The  video  records  were 
analyzed  for  lateral  (spanwise)  streak  spacing,  pocket  width,  and  time  between 
pocket  events,  with  the  streak  spacing  and  pocket  width  normalized  by  w^/v  and 
time  normalized  by  u//v.  Their  results  support  the  proposition  by  Smith  and 
Metzler  (1983)  that  the  normalized  spanwise  streak  spacing  is  a  universal 
number  =  100  and  that  the  spanwise  pocket  scale  (y*  =  yu,/i')  is  related  to  R, 
by 
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Figure  8.  Vertical  profiles  of  inner-tionnalized  streamwise  velocity  statis¬ 
tics:  (a)  mean  velocity,  with  Coles  law  (solid  line)  and  upper  and 

lower  limits  suggested  by  Blackwelder  and  Haritonidis  (1983)  (dotted 
lines);  (b)  velocity  rms,  with  open  symbols  denoting  results  from 
Klewicki  and  Falco  (1990);  (c)  skewness,  with  open  symbols  denoting 
results  from  Klewicki  (1989);  (d)  kurtosis,  with  open  symbols  denoting 
results  from  Klewicki  (1989)  (from  Metzger  and  Klewicki,  1996). 


y*  =  A  log  +  B  , 


(4-2) 


where  A  and  B  are  constants. 

Metzger  and  Klewicki  (1996)  obtained  measurements  of  axial  and  vertical 
velocities  and  velocity  gradients  within  the  buffer  region  and  the  lower 
portion  of  the  inertial  sublayer  using  a  six -element  x-array  probe.  This 
probe  was  positioned  successively  at  0.5,  1.0,  1.5,  and  2.0  m  AGL,  with  10-min 
time  series  collected  at  each  level.  The  fast  sampling  rate  provided  suffi¬ 
cient  data  to  calculate  the  first  four  statistical  moments  of  axial  and 
vertical  velocity  and  their  gradients.  These  measurements  revealed  several 
wall  region  characteristics:  (1)  an  increase  in  velocity  variances  with 
increasing  proximity  to  the  surface,  a  result  consistent  with  the  Falco  (1991) 
and  Willmarth  and  Lu  (1972)  models  which  describe  increased  q*  generation  near 
the  surface;  (2)  an  axial  velocity  skewness  increase  with  proximity  to  the 
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surface,  indicating  a  near-surface  increase  in  intermittent  large  excursions 
from  the  mean;  (3)  a  decrease  in  vertical  velocity  skewness  with  proximity  to 
the  surface,  reaching  a  value  near  zero  at  0.5  m;  (4)  a  near-surface  increase 
in  axial  and  vertical  velocity  gradients  with  proximity  to  the  surface,  with 
even  greater  increases  in  the  kurtosis  of  the  velocity  gradient;  and  (5)  a 
maximum  at  0.5  m,  suggesting  a  buffer  sublayer  origin  for  surface¬ 
generated  Re3molds  stress.  The  increased  near-surface  vertical  velocity 
skewness  and  kurtosis  are  consistent  with  the  presence  of  weak  downward 
"sweeps"  spreading  into  elongated  shear  layers  as  downward  motion  is  blocked 
and  interrupted  by  a  smaller  number  of  energetic  upward-moving  fluid  "ejec¬ 
tions."  Also,  the  vertical  flux  of  Re3molds  stress  is  greatest  at  the  lowest 
measurement  level,  indicating  a  mean  transport  of  turbulent  energy  away  from 
the  surface.  These  preliminary  results  suggest  to  Klewicki  et  al.  (1995)  that 
normalized  stream-wise  turbulence  intensity  and  Re3molds  stress  profiles  may 
fall  on  universal  curves  depicting  the  characteristics  of  near-wall  motions. 
Further  investigations  will  be  required  to  verify  this  hypothesis. 

In  another  SLTEST  site  experiment,  Folz  (1997)  deployed  arrays  of  up  to  24 
P3-rallel  straight  hot-wire  sensors  and  a  modular  12 -sensor  probe  designed  for 
direct  velocity  gradient  tensor  and  shear  measurements  within  the  viscous 
sublayer.  His  major  findings  include  the  stability  and  R,  dependence  of 
several  viscous  sublayer  flow  and  turbulence  characteristics:  (1)  the 
turbulence  intensity  and,  to  a  lesser  extent  skewness,  increase  with  Rj  and 
instability,  but  kurtosis  is  relatively  unaffected  by  these  variables;  (2)  the 
mean  velocity  profile  is  independent  of  R,,  but  the  logarithmic  profile  slope 
is  stability- dependent;  (3)  the  rms  streamwise  and  vertical  velocity  gradients 
are  independent  of  R,,  but  the  gradient  kurtosis  increases  with  R,;  (4) 
turbulent  kinetic  energy  production  and  dissipation  rates  and  vorticity 
production  are  independent  of  R^;  (5)  turbulent  kinetic  energy  q^  dissipation 
is  highly  correlated  with  strong  positive  vertical  and  axial  velocity  gradi¬ 
ents  at  high  R, ,  but  the  vertical  velocity  contributes  little  to  q*  at  low  R,. 

4.4  MEASUREMENTS  AROUND  ISOLATED  CUBES 

Metzger  and  Klewicki  (1996)  provide  a  detailed  description  of  a  series  of 
measurements  made  around  a  0.25-m  cube  positioned  on  the  salt  flat  surface  at 
the  SLTEST  site.  They  also  discuss  measurements  made  in  the  wake  of  the  2-m 
cube  while  it  was  positioned  on  the  ground  near  Tower  Grid  for  the  recircula¬ 
tion  region  studies  (refer  to  Section  2.2  for  details  of  the  Recirculation 
Subtest  trials).  The  purpose  of  the  SLTEST  0.25-m  cube  experiment  was  to 
explore  whether  the  flow  field  around  a  surface -mounted  cube  is  independent  of 
R«,  while  the  2-m  cube  wake  measurements  were  made  to  determine  wake  width  for 
the  Recirculation  Subtest  trials.  Figure  9  shows  the  coordinate  system 
adopted  at  both  sites.  The  origin  is  at  the  center  of  the  bottom  cube  face 
with  X  representing  downwind  distance  (into  the  wake),  y  representing  spanwise 
(crosswind)  distance,  and  z  representing  vertical  distance.  All  distances 
were  normalized  by  cube  face  length  h,  the  characteristic  dimension  of  the 
cube . 

The  SLTEST  cube  trials  included  axial  velocity  and  surface  pressure 
measurements  designed  to  characterize  surface  pressure  and  axial  velocity 
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statistics  in  the  proximity  of  the  0.25-m  cube.  Several  sets  of  pressure  time 
series  measurements  were  taken  at  normalized  spanwise  distances  ranging  from 
0.75  to  2.0.  The  results  show  a  turbulent  kinetic  energy  gain  of  an  order  of 
magnitude  or  more  in  the  disturbed  flow  region,  increasing  with  proximity  to 
the  side  of  the  cube.  Metzger  and  Klewicki  (1996)  conclude  that  large  flow 
obstacles  such  as  the  cube  enhance  the  high  frequency  content  of  the  turbulent 
boundary  layer  through  the  stretching  and  reorientation  of  vorticity.  They 
also  note  that  the  skewness  of  the  pressure  fluctuations  very  near  the  cube 
(y/h  -  0.75,  or  6.25  cm  spanwise  from  the  cube  surface)  is  positive  due  to 
relatively  infrequent  energetic  high  pressure  excursions,  but  becomes  slightly 
negative  at  y/h  >  1  and  remains  so  in  the  free  atmosphere  beyond  this  dis¬ 
tance  . 

Axial  velocity  measurements  were  also  made  at  five  positions  defined  by 
z/h  -  0.125  and  x/h  -  0,  with  y/h  ranging  from  0.6  to  1.24.  Metzger  and 
Klewicki  (1996)  report  that  the  axial  velocity  fluctuation  probability 
distribution  functions  (PDFs)  at  y/h  <1.0  feature  positive  skewness  and  sharp 
peaks,  while  at  y/h  a  1.0  the  flow  begins  to  relax  back  to  that  of  the  undis¬ 
turbed  boundary  layer  with  a  nearly  Gaussian  axial  velocity  PDF.  Test  results 
suggest  that  the  shape,  stability,  and  structure  of  the  horseshoe  vortex  which 
wraps  around  the  cube  is  Re3molds  number  dependent.  Additional  measurements 
are  required  to  verify  this  h3rpothesis. 

Half -wake  widths  behind  the  2-m  cube  were  determined  from  measurements  by 
six  single-element  hot-wire  probes  positioned  at  z/h  -  0.5  and  x/h  -  1.25,  as 
shown  in  Figure  9.  The  purpose  of  these  measurements  was  to  define  the  wake 
width,  which  is  defined  by  Hunt  and  Castro  (1984)  as  the  distance  between 
shear  layer  centers.  The  velocity  profile  shear  layer  center  is  located 
midway  between  the  maximum  and  minimum  velocity  positions  where  the  rms  axial 
velocity  peaks.  A  total  of  eight  successful  10-min  runs  were  accomplished 
consisting  of  six  with  the  front  face  of  the  cube  oriented  normal  to  the  wind 
(face-on)  and  two  with  the  cube  oriented  45®  relative  to  the  mean  wind  (cor¬ 
ner-on)  .  Because  fine  dust  particles  from  the  clay  surface  interfered  with 
the  hot-wire  measurements  in  the  Tower  Grid  area,  fewer  wake  width  measure¬ 
ments  are  available  than  were  planned. 

The  hot-wire  anemometer  measurements  in  the  wake  of  the  2-m  cube  were 
analyzed  for  mean  and  rms  axial  velocity.  Figure  10  shows  the  spanwise  mean 
and  rms  axial  velocity  profiles  for  the  six  face-on  and  two  side-on  cube 
orientation  trials.  Metzger  and  Klewicki  (1996)  use  normalized  distance  to 
the  peak  rms  axial  velocity  from  run  #6  and  assvune  a  symmetric  flow  to 
estimate  a  cube  wake  width  of  y/h  -  1.3  ±0.1  for  the  face-on  orientation, 
which  is  consistent  with  results  from  earlier  studies  at  low  Reynolds  numbers. 
The  side-on  results  (Figures  10c,  lOd)  were  insufficient  to  draw  definitive 
conclusions . 
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Figure  9.  Plan  view  diagram  of  the  2-m  cube  in  its  face-on  (solid  line)  and 
side-on  (dashed  line)  positions  with  respect  to  the  mean  wind,  and 
hot-wire  probe  positions.  Hot-wires  are  positioned  at  a  height 
of  1  m  AGL,  a  distance  of  5h/4  downwind  of  the  cube  center, 
and  spanwise  distances  y/h  of  0.5,  0.75,  1.0,  1.3,  1.7,  and  2.0 
(after  Metzger  and  Klewicki,  1996). 
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SECTION  5.  PLUME  PROFILES  AND  TUBRTIT.F.WrF;  EFFECTS 


5.1  NEAR  SURFACE  PLUME  PROFILE  CHARACTERISTICS 

Vertical  profiles  of  plume  concentrations  were  obtained  during  the  Plume 
Profile  and  Turbulence  Effects  Subtest  conducted  from  21  July  (Julian  date  JJJ 
-=  202)  to  3  August  1995  (JJJ  =  215)  at  the  Tower  Grid  test  site.  Figure  11 
illustrates  the  experimental  configuration.  Each  trial  began  with  the  release 
of  a  known  quantity  of  a  tracer  gas  (propylene)  at  distances  ranging  from  12.5 
to  100  m  upwind  of  a  tower  decorated  with  TIP-SJ2  photoionization  detectors 
(PIDs).  The  gas  was  released  through  a  5-cm  diameter  PVC  pipe,  producing  a 
continuous  propylene  plvune  for  trial  durations  ranging  from  16  to  35  min.  Gas 
release  heights  were  at  0.05,  1.0,  or  2.5  m,  respectively  representing 
releases  in  the  near-wall,  lower,  and  mid  regions  of  the  inertial  sublayer.  A 
triangular  open-mesh  galvanized  metal  9-m  tower  with  a  dimension  of  46  cm  on 
each  side  supported  the  detectors.  The  PIDs  were  mounted  at  heights  ranging 
from  0.25  to  9.0  m,  with  0.5-m  spacing  from  0.5  to  3.0  m  and  1-m  spacing  from 
3.0  to  9.0  m.  Trial  summaries,  with  trial  names  defined  by  the  two-digit 
Julian  date  (JJ)  and  start  times  (nnnn  in  Mountain  Daylight  Time)  are  given  in 
Table  3.  Detailed  plume  profile  information  is  available  in  Appendix  A. 

TIPSJ2  PIDs  have  sufficiently  fine  temporal  and  spatial  response  to 
resolve  the  details  of  a  plume's  internal  concentration  time  history.  Each 
PID  contains  an  orifice ,  an  illumination  chamber  with  an  ultraviolet  lamp ,  a 
set  of  electrically  biased  collector  plates,  a  fan  which  draws  air  through  the 
detector,  and  supporting  electronics.  Air  drawn  through  the  detector  is 
illuminated  by  the  ultraviolet  lamp.  Gases  passing  through  the  illumination 
chamber  that  have  ionization  potentials  of  10.6  electron-volts  (eV)  or  less 
(propylene's  ionization  potential  is  9.73  eV)  are  ionized,  producing  charged 
particles.  Electric  currents  generated  as  the  charged  particles  impinge  upon 
plates  biased  with  the  opposite  charge  produce  voltages.  These  voltages  are 
related  through  calibration  curves  to  airborne  gas  concentrations.  The  TIPSJ2 
frequency  response  is  about  0.0037  s  (6  dB  fall-off  at  270  Hz).  The  TIP-SJ2 
output  was  logged  at  4000  Hz  using  a  fast  sample -and-hold  analog- to -digital 
(A/D)  converter  with  16-bit  resolution.  Further  details  on  the  TIP-SJ2s,  the 
data  collection  system,  and  procedures  used  to  generate  propylene  concentra¬ 
tion  statistics  can  be  found  in  Yee  et  al.  (1993). 

The  plume  profile  summaries  in  Table  3  contain  the  pliime  release  height  Zg, 
plume  centroid  height  z^,  height  of  the  pltune  profile  peak  concentration  Zp, 
peak- to -mean  concentration  ratio  P/M,  and  unconditional  concentration  fluctua¬ 
tion  intensity  i.  The  plume  centroid  height  is  defined  as  the  height  of  the 
detector  with  the  maximum  mean  concentration.  The  peak  concentration  Cp^  at 
each  detector  is  the  concentration  exceeded  1  percent  of  the  time,  and  P/M  is 
the  ratio  of  the  profile -maximum  Cp.  to  the  mean  concentration  at  the  plume 
centroid.  The  unconditional  fluctuation  intensity  i  is  the  standard  deviation 
of  the  measured  concentrations  (zero  concentrations  included)  at  the  plume 
centroid  divided  by  the  mean  plume  centroid  concentration. 
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The  propylene  disseminator  and  tower-mounted  sampling  array  (not 
to  scale) . 


An  understanding  of  the  physical  mechanisms  which  cause  changes  in  plume 
centroid  height  is  important  for  modeling  near-surface  dispersion  because 
model  predictions  are  based  on  the  premise  that  the  plume  centroid  location  is 
known.  In  particular,  Ze  and  Zp  are  assumed  to  be  consistent  with  z,,  with 
changes  in  centroid  height  accounted  for  by  reflection  arguments.  Differences 
between  z, ,  Zp,  and  z^  attributable  to  physical  mechanisms  other  than  surface 
reflection  are  apparent  in  the  Table  3  data.  In  qualitative  terms,  the 
centroids  of  plumes  released  above  the  surface  into  an  adiabatic  atmosphere 
tend  to  descend  toward  the  surface,  while  the  centroids  of  plumes  from  surface 
releases  tend  to  rise.  Also,  Zp  is  typically  found  above  z^.  The  physical 
mechanisms  affecting  the  plume  centroid  height  include  initial  conditions  such 
as  plume  momentum  and  buoyancy  and  "far  field"  effects  such  as  flux  divergence 
and  turbulence  skewness. 
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Biltoft  (1995)  analyzed  the  initial  conditions  for  propylene  released  from 
the  dissemination  system  used  in  these  trials.  He  concluded  that  the  momentrjn 
effect  is  negligible,  but  that  density- driven  slumping  occurs  during  light 
winds.  The  plume  Richardson  number,  defined  by  Chatwin  (1985)  as 

=  g'b/ul  (5-1) 


where  g'  is  a  density-weighted  gravitational  acceleration  and  b  is  a  charac¬ 
teristic  cloud  length,  determines  whether  or  not  material  released  into  the 
atmosphere  will  slump  toward  the  surface.  Chatwin' s  criterion  for  a  gravity- 
dominated  dense  gas  is  R,,  »  0.2;  neutrally -buoyant  behavior  is  expected  for 
Rp  «  0.2,  and  a  transitional  condition  exists  for  Rp  on  the  order  of  0.2. 
Using  these  criteria,  propylene  typically  exits  the  disseminator  as  a  dense 
gas,  but  entrainment  of  ambient  air  into  the  plume  rapidly  reduces  Rp  below 
the  slumping  threshold  unless  the  release  occurs  during  quiescent  nocturnal 
conditions  (u.  s  0.1  m/s).  Under  quiescent  conditions,  denser-than-air  gases 
such  as  propylene  slump  toward  the  surface  and  fan  out  in  a  shallow  pool  until 
sufficient  mixing  occurs  to  overcome  gravitational  effects.  Slumping  likely 
lowered  the  centroid  heights  for  Trials  P030545,  P100550,  P100705,  P122325, 
P132250,  P132340,  P140030,  P142154,  and  P142245  of  this  series. 

In  addition  to  the  initial  conditions,  the  plume  centroid  position  is 
determined  by  characteristics  of  the  turbulent  wind  field  into  which  it  is 
released.  Raupach  (1983)  considers  the  dispersion  of  a  neutrally -buoyant 
tracer  into  an  idealized  adiabatic  turbulent  flow  and  derives  an  analytic 
expression  for  the  mean  tracer  particle  vertical  drift  velocity  u  as  a 
function  of  u..  This  equation  is 

u(t)  =  jbu,  (1 1  +  t) )  ,  (5-2) 


where  b  is  a  constant  on  the  order  of  0.4  and  t  is  travel  time  non- dimens iona- 
lized  by  the  Lagrangian  time  scale.  Thus,  the  vertical  plume  centroid 
position  is  related  to  measureable  ABL  quantities  such  as  u,  and,  through 
Equation  (3-12),  to  Superimposed  on  this  idealized  relationship  are  non- 

Gaussian  effects  such  as  skewed  turbulence,  which  enhances  plume  centroid  rise 
during  convection  (see  Wyngaard  and  Weil,  1991),  and  stable  thermal 
stratification,  which  suppresses  plume  centroid  rise. 

The  plume  centroid  height  also  influences  the  magnitude  of  the  concentra¬ 
tion  fluctuation  intensity  i  and  peak-to-mean  ratio  P/M.  The  fluctuation 
intensities  presented  in  Table  3  are  those  found  at  the  plume  centroid.  These 
are  "total"  or  "unconditional"  statistics,  which  include  the  zero  concentra¬ 
tions  readings  that  occurred  when  the  plume  meandered  away  from  the  detector 
array.  Yee  et  al.  (1993)  suggest  that  i  at  the  plume  centroid  can  be  repre¬ 
sented  as  a  power  law  function  of  downwind  travel  distance  normalized  by  z^. 
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i  =26.6  (x/z^) 


(5-3) 


Implicit  in  Equation  (5-3)  is  the  assumption  that  the  plume's  height  above  the 
surface  largely  determines  the  intensity  of  small-scale  turbulence  acting  on 
the  plume.  For  example,  a  plume  released  near  the  surface  will  be  subjected 
to  more  intense  small  scale  mixing  (reducing  i)  than  a  plume  released  above 
the  surface.  However,  measured  fluctuation  intensities  exhibit  considerable 
scatter  about  the  fluctuation  intensities  calculated  using  Equation  (5-3). 

The  fluctuation  intensity  is  often  not  well  defined  using  the  x/z^  normaliza¬ 
tion  because  the  centroid  of  a  plume  released  into  the  atmosphere  does  not 
necessarily  remain  at  the  release  height,  and  the  plume  is  subjected  to 
decreased  (or  increased)  internal  mixing  depending  on  whether  the  centroid 
rises  (or  sinks).  A  considerable  decrease  in  scatter  was  obtained  with  Equa¬ 
tion  (5-3)  using  an  x/z^  normalization,  suggesting  that  the  plume  centroid 
vertical  position  time  history  must  be  known  to  accurately  estimate  i. 

The  assumption  that  the  plume  peak  concentration  is  found  at  the  plume 
centroid  is  frequently  violated,  particularly  when  Zc  originates  near  the 
surface.  The  intensity  of  fine  scale  turbulence  is  at  a  maximum  near  the 
surface  (see  Stull,  1988).  Thus,  the  concentration  gradients  within  the  near¬ 
surface  portion  of  the  plume  are  rapidly  destroyed  by  intense  mixing.  This 
mixing  reduces  both  the  concentration  peaks  and  adjacent  pockets  of  clear  air, 
while  leaving  the  mean  concentration  relatively  unchanged.  Near  surface 
mixing  operates  on  a  timescale  shorter  than  the  timescale  of  the  turbulent 
fluxes  that  elevate  the  plume  centroid.  The  result  is  that  Zp  often  remains 
at  or  rises  above  Zg,  particularly  when  the  plume  centroid  sinks  toward  the 
surface . 

In  spite  of  its  importance  to  dispersion  modeling,  only  a  rudimentary 
understanding  of  near-surface  plume  profile  behavior  exists.  While  studies 
like  those  performed  by  Raupach  (1983)  and  Wyngaard  and  Weil  (1991)  provide 
useful  insights  into  pltjune  centroid  behavior,  the  development  of  a  physically 
valid  plume  centroid  vertical  displacement  model  will  require  improved 
understanding  of  scalar  variance  vertical  transport  and  its  relationship  to 
quantities  such  as  the  fluxes  of  heat  and  momentum.  Specifically,  the 
contributions  of  the  velocity- concentration  correlation  and  correlation 
gradient  terms  in  the  scalar  variance  transport  equation  must  be  better 
understood,  which  will  require  the  development  of  instrumentation  capable  of 
providing  velocity- concentration  correlation  measurements. 

5.2  NEAR  SURFACE  WIND  AND  TURBULENCE  ANALYSES 

The  vertical  array  of  PIDs  at  the  Tower  Grid  test  site  was  accompanied  on 
the  9-m  tower  by  Applied  Technologies,  Inc.  sonic  anemometer/ thermometers 
(sonics),  each  consisting  of  an  ultrasonic  sensor  array  and  an  electronics 
module.  The  sonic  arrays  were  mounted  with  an  orientation  toward  240°  with 
respect  to  true  North,  Two-axis  sonic  arrays  (Model  RSWS-201/2A)  were  located 
at  the  1.5-  and  6.0-m  levels,  and  tri-axis  sonics  (Model  RSWS-201/3A)  were 
located  at  the  3.0-  and  9.0-m  levels.  The  two-axis  sonics  provided  measure¬ 
ments  of  the  horizontal  wind  components  which  were  rotated  into  the  mean 
alongwind  u  and  crosswind  v  components  during  data  processing.  The  tri-axis 
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sonics  provided  measurements  of  all  three  components  of  the  wind  vector  (u,  v, 
and  vertical  w)  which  were  subjected  to  an  identical  horizontal  coordinate 
rotation  into  the  mean  wind.  All  of  the  sonics  also  measured  the  speed  of 
sound,  which  was  subsequently  converted  into  sonic  temperature  (nearly 
identical  to  virtual  temperature) .  The  sonic  data  collected  at  a  rate  of  10 
Hz  provided  sufficient  temporal  and  spatial  resolution  to  characterize  the 
surface  boundary  layer  turbulence  in  the  vicinity  of  the  tower.  Statistical 
summaries  of  the  sonic  data  collected  during  dissemination  periods  are 
presented  in  Appendix  A. 

Disseminations  for  the  Plume  Profile  and  Turbulence  Effects  Subtest  were 
planned  for  wind  speeds  ranging  from  1  to  6  m/s  and  wind  directions  within 
±60®  of  180°.  Most  of  the  trials  occurred  with  wind  speed  less  than  3  m/s. 

An  attempt  was  made  to  select  dissemination  time  periods  when  the  wind 
direction  was  steady,  but  abrupt  changes  in  wind  direction  due  to  turbulence 
bursts  or  terrain- induced  meander  often  carried  the  plume  away  from  the  array. 
Attempts  were  made  during  some  trials  to  realign  the  source  to  increase  plvune 
passage  through  the  detector  array,  but  most  data  sets  contain  large  blocks  of 
zero  concentrations  due  to  persistent  meander. 

Table  4  gives  the  Plume  Profile  and  Turbulence  Effects  Subtest  trial 
micrometeorological  summaries.  These  summaries  include  wind  and  temperature 
statistics,  friction  velocity,  and  the  Obukhov  length  as  a  stability  para¬ 
meterization.  Because  none  of  the  eddy  correlation- derived  friction  veloci¬ 
ties  met  the  Busch  and  Panofsky  (1968)  quality  criterion  (u.  >0.32  m/s). 
Equation  (3-12)  was  used  to  estimate  u..  Table  4  also  provides  uV  computed 
from  sonic  measurements  at  the  1.5-  and  6-m  levels,  and  the  Diabatic  Ratio, 
which  was  derived  from  covariance  quadrant  contributions  to  uV  (discussed 
below) . 

An  interesting  observation  from  this  trials  series  is  that,  although 
detailed  meteorological  measurements  were  made  near  the  times  of  transition 
through  sunrise  and  sunset,  conventional  "neutral"  stability  was  not  observed. 
Both  the  Obukhov  length  and  Diabatic  Ratio  (defined  later)  indicate  that, 
during  the  transition  from  a  very  unstable  to  a  very  stable  ABL,  both  stabili¬ 
ties  exist  simultaneously  within  different  layers.  For  example,  measurements 
from  the  1.5-  and  6.0-m  levels  during  Trial  P122325  (Table  4)  indicate  a  very 
stable  layer  within  a  few  meters  of  the  surface  with  an  unstable  layer 
persisting  above  it  well  into  the  night.  It  appears  that  the  nocturnal  stable 
layer  develops  near  the  surface  and  gradually  progresses  through  successively 
deeper  ABL  layers  after  sunset.  Thus,  transition  in  the  desert  ABL  consists 
of  a  superposition  of  diabatic  layers  rather  than  a  true  adiabatic  ABL.  The 
extreme  stability  conditions  indicated  by  the  small  Obukhov  lengths  calculated 
for  many  of  the  trials  are  artifacts  of  very  light  winds,  which  produce  near- 
zero  (and  statistically  unstable)  fluxes. 

Another  characteristic  of  the  Table  4  trial  summaries  is  the  extreme 
stability  conditions ,  as  indicated  by  the  small  Obukhov  lengths  calculated  for 
many  of  the  trials .  These  small  Obukhov  lengths  are  artifacts  of  the  light 
winds  which  produce  small  values  of  and  u.  which,  when  cubed,  produce  a 
small  L.  The  abrupt  change  between  similarity-based  stability  extremes  that 
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No  data  available. 


occurs  during  the  high  desert  boundary  layer  transition  through  sunrise  and 
sunset  in  light  winds  is  not  accompanied  by  comparable  changes  in  turbulence 
levels,  and  hence  by  extreme  differences  in  the  rate  of  diffusion.  This 
situation  calls  into  question  the  adequacy  of  similarity-based  boundary  layer 
representations,  particularly  at  non- ideal  sites  and  during  non- steady 
conditions  where  adequate  vertical  velocity  measurements  or  reliable  u. 
estimates  are  unlikely  to  be  available. 

The  limitations  of  similarity-based  boundary  layer  representations  suggest 
a  need  for  a  more  robust  method  for  describing  the  turbulent  state  of  the 
atmosphere,  particularly  during  non- ideal  conditions.  Desirable  characteris¬ 
tics  for  an  improved  stability  and  turbulence  parameter  include:  (1)  ease  of 
calculation  in  near  real-time  from  near-surface  measurements  using  a  single 
robust  instrument,  (2)  sensitivity  to  changes  in  stability  and  turbulence 
levels,  (3)  statistical  stability  during  non- ideal  conditions,  and  (4)  a 
mathematically  tractible  and  physically  meaningful  scale  range.  In  this 
study,  the  search  for  a  new  stability  parameterization  focused  on  the  IFF 
covariance  term  in  the  stress  budget  equation  (see  term  [3]  of  Equation  (3-1)) 
because  Indices  based  on  velocity- temperature  covariance  measurements  from  a 
two-axis  sonic  anemometer/thermometer  appear  to  most  nearly  satisfy  these 
criteria. 


The  u't'  covariance  offers  several  practical  advantages  over  the  and 
IFF  covariances  as  a  turbulence  indicator:  (1)  the  absence  of  a  vertical 
velocity  measurement  requirement  relieves  the  need  for  research -grade  instru¬ 
ment  (such  as  a  tri-axis  sonic)  and  precise  sensor  leveling  and  alignment  (for 
a  discussion  of  these  effects,  see  Haugen  and  Kaimal,  1969);  (2)  measurements 
made  only  in  the  horizontal  plane  can  be  obtained  from  two-axis  sonic  arrays 
positioned  close  (within  a  meter,  depending  on  Ze)  to  the  surface,  while  tri- 
axis  sonics  should  be  mounted  at  higher  levels;  (3)  iFF  is  statistically  more 
stable  than  u'w'  (a  10-min  sampling  period  appears  to  provide  a  reasonable  u'T' 
estimate)  ;  and  (4)  because  u'  is  usually  larger  than  w'  ,  uV  is  typically 
larger  than  the  uV  or  (FF  obtained  during  the  same  sampling  period. 


The  PF  covariance  is  t3rpically  considered  a  horizontal  temperature  flux 
term  (see,  for  example,  Wyngaard  et  al.,  1971),  but  its  magnitude  and  sign 
appear  to  vary  with  stability  and  with  day-night  transitions  in  ways  that  are 
similar  to  the  variations  in  iFF.  Also,  except  in  the  presence  of  dynamic 
frontal  boundaries,  the  horizontal  temperature  gradient  does  not  support  the 
observed  magnitude  of  PF  (i.e.,  horizontal  temperature  advection  is  usually 
weak) .  Because  vertical  momentum  and  temperature  gradients  usually  exceed 
their  horizontal  gradients  by  several  orders  of  magnitude,  the  fluctuating 
velocity  and  temperature  components  that  make  the  predominant  contribution 
to  PF  most  likely  originate  with  vertical  eddy  exchanges.  That  is,  the 
vertical  flux  of  the  temperature -momentvun  covariance  is  likely  to  contribute 
more  than  the  horizontal  temperature  flux  to  the  magnitude  of  PF.  The 
remainder  of  this  section  examines  the  relationship  of  PF  to  non- steady 
boundary  layer  processes  using  Fourier  spectrum  analysis  and  covariance 
quadrant  analysis. 
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Fourier  spectrum  analysis  involves  the  conversion  of  time  series  data  from 
the  time  domain  into  the  frequency  domain.  Fourier  analysis  is  based  on  the 
premise  that  any  continuous  time  series  can  be  represented  by  a  finite  number 
of  sinusoidal  components  sximmed  over  a  finite  frequency  range.  It  offers  the 
advantage  of  describing  the  distribution  of  the  variance  or  covariance  of 
time -varying  quantities  among  the  range  of  frequency  scales  that  form  the  time 
series.  That  is,  Fourier  analysis  partitions  energy  by  frequency  rather  than 
by  time.  The  Fourier  transformation  partitions  power  (change  in  a  variable 
per  unit  time)  of  a  time  series  into  harmonic  frequency  components.  An 
advantage  offered  by  this  power  spectrum  partition  is  that  the  power  (or 
turbulent  energy)  within  one  harmonic  is  independent  of  the  amplitudes, 
phases,  and  frequencies  of  other  harmonics  present  within  the  time  series. 
Thus,  spectrum  analysis  displays  the  distribution  of  turbulent  energy  content 
within  a  time  series  as  a  function  of  turbulence  scale.  Additional  informa¬ 
tion  on  time  series  analysis  for  boundary  layer  phenomena  is  given  by  Stull 
(1988)  and  Kaimal  and  Finnigan  (1994). 

Relationships  between  two  time  series  can  be  evaluated  by  examination  of 
the  complex  products  of  their  spectra  such  as  the  cross  spectrum,  phase,  and 
coherence.  The  cross  spectrum  includes  the  cospectrum  (the  in-  or  out-of- 
phase  portion  of  two  spectral  frequencies)  and  the  quadrature  (that  portion  of 
the  paired  frequencies  shifted  by  Va  wavelength,  or  90®).  The  phase  and 
coherence  provide  the  same  information  in  a  normalized  polar  coordinate 
system,  with  the  angular  phase  relationship  defined  by 

Phase  ■=  ATAN2 (quadrature/cospectra)  (5-4) 

where  ATAN2  is  an  arctangent  function  expanded  to  the  range  -tt  to  tt.  Phase 
information  is  presented  in  degrees ,  with  0°  indicating  an  in-phase  and  180° 
representing  an  out-of-phase  relationship.  The  coefficient  of  coherence  is  a 
measure  of  the  correlation  between  two  time  series  presented  as  a  function  of 
frequency.  It  is  given  by 

Coherence  =  (COS^  +  QUAD^)/(Spctrl)*(Spctr2) ,  (5-5) 

where  Spctrl  and  Spctr2  are  the  spectral  component  magnitudes  of  time  series  1 
and  2  and  COS  and  QUAD  respectively  denote  the  cospectrum  and  quadrature 
spectrum.  Kaimal  and  Gaynor  (1983)  describe  the  fast  Fourier  transform  (FFT) 
spectrum  analysis  program  used  to  process  the  sonic  10-Hz  u  and  T  data.  The 
spectra  produced  by  this  FFT  program  are  smoothed  through  block  averaging, 
tapered,  scaled  to  meter -kilogram- second  (mks)  units,  and  normalized  by 
multiplying  each  harmonic  component  by  its  frequency.  Fourier  analysis 
of  u'T'  provides  information  on  how  the  surface  layer  responds  to  the  influx  of 
energy  during  a  transition  through  sunrise,  as  shown  in  the  Figure  12  and 
summarized  in  Table  5.  The  details  of  Figure  12  and  Table  5  are  discussed  as 
part  of  the  case  study  analysis  in  Section  5.3. 
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Another  technique  applicable  to  uV  studies  is  covariance  quadrant  analy¬ 
sis.  The  covariance  of  any  quantities  A  and  B,  where  A'  and  B'  represent 
fluctuations  about  the  means  A  and  B,  include  positive  (greater  than  the 
mean)  and  negative  (smaller  than  the  mean)  A'  and  B'  values.  Thus,  there  are 
four  possible  combinations  of  A'  and  B' :  (1)  A'>0,  B'>0;  (2)  A'>0,  B'<0;  (3) 
A'<0,  B'<0;  and  (4)  A'<0,  B'>0.  If  either  A'  or  B'  is  near-zero,  the  contri¬ 
bution  of  the  A'B'  product  to  the  flux  is  near-zero.  Thus,  the  covariance  can 
be  divided  into  four  quadrants  which  contribute  to  the  flux  plus  a  fifth  no- 
flux  region  near  the  zero  axes  (see  Figure  13) .  Because  the  velocity  resolu¬ 
tion  threshold  for  ATI  sonics  is  near  0.03  m/s,  this  threshold  was  defined  as 
the  boundary  of  the  no-flux  region.  As  noted  by  Raupach  (1981),  the  stun  of 
the  covariance  fraction  contributions  from  all  of  the  quadrants  is  the  total 
covariance.  Covariance  quadrants  and  the  no-flux  region  for  IPf'  are  depicted 
in  Figure  13. 


The  similarities  in  the  diurnal  variations  of  'iFt'  and  w'f’  can  be  explained 
by  a  simple  model  of  the  boundary  layer.  Consider  a  horizontally  homogeneous 
boundary  layer  in  which  discrete  eddies  move  through  a  horizontal  plane  in 
response  to  a  Gaussian  vertical  velocity  distribution  with  zero  mean.  Because 
wind  speed  increases  with  height  in  this  idealized  boundary  layer,  eddies 
descending  from  above  the  horizontal  plane  produce  wind  speed  fluctuations  u' 
above  the  mean  wind  speed  in  the  plane,  while  eddies  rising  from  below  produce 
wind  speed  fluctuations  below  the  mean  wind  speed.  Thus,  u'-H  is  associated 
with  a  negative  vertical  velocity  fluctuation  w' -  and  u' -  is  associated  with  a 
positive  vertical  velocity  fluctuation  w'+.  Quadrant  Q1  (u'+,  T'+)  and 
quadrant  Q3  (u'-,T'-)  eddies  moving  through  the  horizontal  plane  both  make 
positive  contributions  to  1P¥.  Following  the  reasoning  given  above,  the  signs 
on  the  horizontal  velocity  fluctuations  indicate  that  the  Q1  eddies  are 
associated  with  the  downward  flux  of  air  that  is  warmer  than  the  air  at  the 
height  of  the  horizontal  plane  and  the  Q3  eddies  are  associated  with  the 
upward  flux  of  cooler  air.  The  stable  thermal  stratification  implied  by  the 
Q1  and  Q3  eddies  suggests  that  they  are  driven  by  shear- induced  turbulence 
rather  than  convection.  Similar  reasoning  suggests  that  the  Q2(u'+,T'-)  and 
Q4(u'-,T'+)  eddies,  which  make  negative  contributions  to  u't',  are  respectively 
associated  with  the  downward  flux  of  cooler  air  and  upward  flux  of  warmer  air. 
The  unstable  stratification  implied  by  the  Q2  and  Q4  eddies  is  consistent  with 
convection. 
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The  simple  boundary  layer  model  discussed  above  leads  to  the  hypothesis 
that  the  magnitudes  of  the  positive  (Quadrants  Q1  and  Q3)  and  negative 
(Quadrants  Q2  and  Q4)  contributions  to  the  total  can  be  used  to  estimate 

the  relative  contributions  of  shear-driven  and  convection- driven  turbulence. 
A  new  stability  parameter,  the  "Diabatic  Ratio,"  is  then  given  by 

DR  =  \Q2  *  Q4\/  {Q1  +  Q3)  (5-6) 


where  Q1 ,  Q2 ,  Q3 ,  and  Q4  are  the  contributions  from  the  four  quadrants  to 
u'T'.  (Note  that  the  numerator  of  Equation  (5-6)  is  the  absolute  value  of  the 
sum  of  the  quadrants  with  negative  contributions  to  uV. )  DR  serves  as  a 
simple  stability  indicator  with  DR<1  indicating  shear -dominated  turbulence  and 
DR>1  indicating  convection- dominated  turbulence.  If  DR*1,  neither  shear  nor 
convection  dominates,  which  is  consistent  with  an  adiabatic  boundary  layer. 

The  simple  boundary  layer  model  also  suggests  that  the  magnitude  of  the 
contribution  in  each  quadrant  is  directly  proportional  to  the  vertical 
velocity  fluctuations  that  move  eddies  through  the  horizontal  measurement 
plane.  Thus,  the  "Total  Diabatic  Influence"  given  by 

TDI  =  \0i\  (5-7) 


should  be  proportional  to  the  mean  vertical  velocity  variance  This 

hjrpothesis  was  independently  tested  using  36  tri-axis  sonic  anemometer  data 
sets .  Part  of  the  data  sets  came  from  sonics  mounted  at  2.0  m  AGL  over  salt 
flats  in  the  Great  Salt  Lake  Desert  (Biltoft,  1997)  and  part  came  from  sonics 
mounted  at  4.0  m  AGL  in  the  San  Joaquin  Valley  in  California.  Wind  speeds  on 
the  salt  flats  were  light  to  moderate  (2.1  to  8.4  m/s),  while  the  San  Joaquin 
Valley  winds  were  very  light  (0.6  to  1.6  m/s).  Stabilities  ranging  from  very 
stable  to  very  unstable  were  represented  in  this  integrated  data  set.  Trial 
durations  varied,  but  were  typically  60  min  or  longer.  No  trend  removal 
filter  was  applied  to  these  data.  The  linear  correlation  between  TDI  and  the 
measured  ^  was  0.835,  indicating  that  70  percent  of  the  vertical  velocity 

variance  was  explained  by  its  relationship  to  TDI.  This  result  is  significant 
based  on  an  F  test  at  the  1-percent  level.  Section  5.3  below  demonstrates  the 
utility  of  DR  and  TDI  for  characterizing  stability  and  boundary  layer  turbu¬ 
lence  during  non- steady  conditions. 

Several  refinements  to  the  DR  and  TDI  are  possible.  The  present  method 
gives  equal  weights  to  the  contributions  from  all  four  quadrants  without 
consideration  of  the  number  of  occurrences  in  each  quadrant.  Weighting  each 
covariance  contribution  by  the  percent  of  occurrences  within  that  quadrant 
could  improve  the  relationship  of  TDI  to  Also,  the  contribution  from  the 

no -flux  region  has  been  ignored  because  its  contribution  to  the  flux  is 
minimal  except  during  very  light  wind  conditions .  A  proportionality  constant 
of  4  (i.e.,  «  TDI/4)  appears  to  apply  in  strongly  diabatic  boundary  layers 
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where  transitions  through  neutral  are  rapid,  but  may  not  work  well  in  situa¬ 
tions  where  the  atmosphere  has  had  sufficient  time  to  become  well  mixed,  as  in 
the  persistent  cloudy  conditions  of  Northern  Europe.  Further  examination  of 
DR  and  TDI  applicability  as  stability  and  turbulence  indicators,  including 
covariance  quadrant  contributions  to  the  Reynolds  equations,  is  beyond  the 
scope  of  this  report,  but  should  be  undertaken  as  follow-on  tasks. 

5.3  CASE  STUDIES  OF  SUNRISE  TRANSITION  AND  NOCTURNAL  TURBULENCE  INTERMITTENCY 

A  series  of  four  trials  on  22  July  1995  were  selected  for  a  case  study 
analysis  of  the  transition  through  sunrise.  This  day  was  selected  because  it 
is  representative  of  a  high  desert  sunrise  transition  from  quiescent  nocturnal 
conditions  with  wind  speeds  on  the  order  of  1  m/s  to  a  period  of  strong 
convection  in  light  winds.  Representative  10-min  periods  were  selected  from 
sequential  trials  to  document  conditions  during  this  transition.  The  10-min 
time  blocks  represent  a  compromise  between  the  need  for  a  sample  that  is 
sufficiently  long  for  second-order  statistics  to  be  statistically  stable,  but 
sufficiently  short  to  be  free  of  trends  and  meander.  The  micrometeorological 
simimaries  for  each  of  these  10-min  trial  segments  are  given  in  Table  6. 

The  micrometeorological  stammaries  (Table  6)  of  wind  and  temperature  mea¬ 
surements  made  within  the  first  2  m  of  the  surface  illustrate  changes  in  the 
wind  and  temperature  fields  during  a  morning  sunrise  transition.  Temperature 
readings  were  obtained  from  the  fiberoptic -quartz  thermometer  with  probes  at 
the  surface  and  2.0  m  AGL.  A  2 -axis  sonic  provided  horizontal  wind  component 
and  sonic  temperature  fluctuation  measurements  at  1.5  m,  while  a  3-axis  sonic 
supplied  vertical  velocity  measurements  at  3.0  m  AGL.  The  mean  wind  speed 
remained  light,  in  the  range  of  1  to  2.5  m/s,  with  the  wind  direction  from  the 
southwest  for  the  entire  transition  period.  The  surface  and  2.0-m  tempera¬ 
tures  remained  fairly  constant  prior  to  sunrise,  and  then  increased  rapidly 
after  0709  MDT  when  direct  sunlight  first  reached  the  ground  surface  on  the 
grid  (sun  on  the  grid  is  delayed  by  the  shadow  of  Camel  Back  Ridge).  The  rate 
of  surface  temperature  increase  initially  lagged  that  of  the  air  immediately 
above  the  surface,  possibly  due  to  heat  absorption  by  the  cool  underlying 
surface  acting  as  a  thermal  sink.  By  the  last  time  period  (0901  MDT) ,  the 
surface  thermal  lag  had  been  overcome  and  the  ground  temperature  had  rapidly 
increased  to  33.5  °C.  It  remained  10  to  15  °C  above  that  of  the  overlying  air 
for  the  remainder  of  the  trial  period. 

The  arrival  of  sunlight  on  the  grid  at  0709  MDT  produced  an  immediate 
increase  in  a^,  a„,  a^,  and  FF,  but  the  DR  for  the  0700-0710  time  period 
remained  at  0.08  due  to  a  relatively  weak  contribution  from  the  second  and 
third  quadrant  covariance  terms.  The  weak  contributions  from  these  quadrants 
indicates  the  absence  of  convection,  even  though  the  heating  produced  enough 
energy  to  increase  the  general  level  of  turbulence  by  a  factor  of  two  to 
three.  The  onset  of  convection,  as  indicated  by  a  sign  change  in  u't'  and  an 
increase  in  the  DR  to  2.22,  had  occurred  by  the  0750  MDT  time  period.  Convec¬ 
tion  became  well  established  between  0750  and  0901,  with  the  DR  contributions 
from  the  second  and  third  quadrants  exceeding  those  of  the  first  and  fourth 
quadrants  by  a  factor  of  two. 
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A  detailed  picture  of  turbulence  development  through  the  sunrise  transi¬ 
tion  is  available  through  analysis  of  the  temperature -velocity  cospectra 
measured  at  1.5  m  AGL,  as  presented  in  Table  5.  Thirty-minute  time  periods 
beginning  at  0645,  0743,  and  0845  MDT  were  selected  for  spectrum  analysis. 

The  30-min  periods  were  selected  to  obtain  statistically  stable  samples.  The 
pre- sunrise  period  beginning  at  0645  MDT  featured  positive  velocity- tempera¬ 
ture  cospectra  with  phase  angles  in  the  range  of  360±60°  (indicating  a  strong 
in-phase  relationship  between  u'  and  T'  over  most  of  the  spectrum).  A  phase 
change,  with  negative  cospectra  and  phase  angles  shifted  to  180+60°  is  shown 
Table  5  for  the  0743  MDT  period.  This  out-of-phase  relationship  between  u' 
and  T'  (consistent  with  the  increase  in  covariance  Quadrants  Q2  and  Q3)  is 
characteristic  of  incipient  convective  layer  development.  Also  evident  from 
the  phase  angle  time  sequence  is  that  convective  layer  development  begins  at 
the  higher  frequencies  and  gradually  expands  into  the  lower  frequencies  with  a 
transition  through  quadrature.  This  transition  is  apparent  in  the  0743  MDT 
data  set  which  contains  a  significant  portion  of  low  frequency  spectral  energy 
in  quadrature  with  phase  angles  of  090±60°.  The  transition  stage  is  followed 
by  full  convection  represented  by  the  0845  MDT  period,  where  phase  angle 
scatter  has  become  more  random,  but  remains  centered  in  the  southeast  quad¬ 
rant. 

Comparable  covariance  and  spectrum  analyses  performed  at  the  3-  and  6-m 
levels  (not  shown)  reveal  the  sunrise  transition  beginning  in  the  buffer 
sublayer  and  progressively  expanding  upward  through  the  boundary  layer.  The 
0743  MDT  time  interval  represents  a  mid- transition  period,  with  the  lowest  3  m 
of  the  boundary  layer  experiencing  the  effects  of  convection  while  the  layers 
above  remain  stable.  The  result  is  a  DR  of  2.22  at  1.5  m  and  0.46  at  6-m 
during  this  time  period,  indicating  a  diabatic  superposition  (stable  over 
unstable) .  This  suggests  that  stable  and  unstable  regimes  coexist  within  a 
few  meters  of  each  other  in  a  diabatic  boundary  layer  undergoing  sunrise 
transition.  In  this  scenario,  a  true  adiabatic  condition  exists  only  momen¬ 
tarily  at  the  level  within  the  ABL  where  the  stability  transition  is  occur¬ 
ring;  an  adiabatic  state  is  not  representative  of  the  entire  ABL  at  any  time 
during  the  transition. 

The  effects  of  sunrise  on  vertical  velocity  spectra  at  3-m  AGL,  as  shown 
in  Figure  12,  provide  additional  insights  into  a  boundary  layer  in  transition. 
The  0545  MDT  spectrum  (A- -A  in  Figure  12)  is  t3rpical  of  nighttime  conditions 
with  a  distinct  low  energy  "spectral  gap"  at  frequencies  below  0.1  Hz.  The 
increased  energy  present  at  sunrise  (0645  MDT,  B--B)  is  observed  across  all 
spectral  bands  and  is  particularly  pronounced  in  the  spectral  gap.  The 
greatest  energy  increase  occurs  within  the  hour  after  the  sun  reaches  the  grid 
(0745,  C--C).  A  familiar  convective  boundary  layer  spectrum  with  a  fully 
developed  energy -containing  region  and  inertial  sublayer  is  present  by  0845 
MDT  (D--D).  An  interesting  feature  of  this  sequence  of  spectrum  plots  is  the 
presence  of  a  high  frequency  energy  "overshoot"  rather  than  an  inertial 
subrange  in  the  0745  MDT  spectra.  One  hypothesis  is  that  this  high  frequency 
energy  overshoot  is  a  consequence  of  the  restriction  on  convective  motions  by 
the  inversion  "lid."  Just  as  a  lid  increases  the  intensity  of  boiling  in  a 
pot,  the  inversion  lid  on  a  developing  convective  layer  prevents  convection 
from  achieving  its  full  vertical  development.  The  result  is  increased  energy 
available  for  convective  mixing  at  the  higher  frequencies  within  the  confines 
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of  the  inversion  lid.  This  confined  band  of  highly  energetic  turbulence 
should  enhance  dispersion  within  the  developing  convective  layer,  producing  a 
well  mixed  plume  within  a  shallow  layer  just  above  the  surface.  Therefore, 
sunrise  transition  can  produce  transient  episodes  of  very  high  concentrations 
near  the  surface  if  the  growing  convective  layer  penetrates  through  stable 
layers  above  the  surface  which  contain  a  trapped  plume. 

In  contrast  to  the  sunrise  transition,  the  turbulence  burst  is  an  inter¬ 
mittent  nocturnal  event  that  occurs  as  a  consequence  of  hydrodynamic  instabil¬ 
ity  within  a  stably- stratified  atmosphere.  Businger  (1973)  offers  the 
following  qualitative  description  of  the  process  in  which  sporadic  turbulence 
events  occur  within  an  otherwise  quiescent  nocturnal  boundary  layer.  When  the 
Flux  Richardson  nxamber  (Equation  (3-14))  exceeds  its  critical  value  (Rp  > 

0.2),  turbulence  is  suppressed  and  laminar  layers  develop.  This  stratifica¬ 
tion  interrupts  the  downward  flux  of  heat  and  momentum.  Near  the  surface,  the 
winds  become  calm  while  winds  aloft  accelerate  because  their  momentvim  is  not 
being  dissipated  toward  the  surface.  The  acceleration  of  wind  accompanied  by 
weak  heat  flux  eventually  drives  Rp  below  its  critical  value,  whereupon  the 
laminar  layers  burst  into  turbulence,  creating  a  sudden  exchange  of  heat  and 
momentiam.  Quiescent  conditions  return  with  the  dissipation  of  this  energy  and 
the  process  repeats  itself. 

Measurements  at  DPG  during  periods  with  light  winds  indicate  that  the 
nocturnal  boundary  layer  in  a  high  desert  environment  is  more  complicated  than 
a  series  of  stratified  layers  intermittently  disturbed  by  turbulence  bursts 
tripped  by  a  Rp  switch.  For  example,  the  20-min  time  series  plots  of  wind 
components  and  speed  of  sound  shown  in  Figure  14  for  Trial  P142154  illustrate 
a  quiescent  nocturnal  light  wind  period  followed  by  a  turbulence  burst.  The 
time  series  show  a  gradual  decay  in  turbulence  level  for  the  first  10  min, 
followed  by  a  second  10 -min  period  that  begins  with  abrupt  perturbations  in 
velocity  components  and  the  speed  of  sound.  Statistical  summaries  for  these 
two  successive  10-min  time  blocks  (identified  by  Julian  date  and  start  times 
2142205  and  2142215)  are  presented  in  Table  7.  The  turbulence  burst  is 
characterized  by  large  increases  in  the  velocity  component  and  temperature 
variances.  This  event  also  caused  the  TDI  to  increase  by  a  factor  of  3,  while 
the  DR  remained  small  due  to  weak  contributions  from  u'T'  Quadrants  Q2  and  Q3. 

A  high  TDI  combined  with  a  low  DR  charcterizes  this  turbulence  burst  as  an 
energetic,  shear -driven  event. 

The  turbulence  burst  that  occurred  during  Trial  P142154  had  a  dramatic 
effect  on  propylene  gas  concentrations  measured  on  the  9-m  tower  sampling 
array.  High  gas  concentrations  were  measured  at  all  tower  levels  between  2205 
and  2215  MDT.  The  concentration  signal  abruptly  disappeared  as  the  turbulence 
event  began,  with  the  subsequent  appearance  of  only  weak,  intermittent 
concentrations  for  the  remainder  of  the  trial.  Concentration  statistics  for 
Trial  P142154  are  based  primarily  on  concentrations  measurements  prior  to  2215 
MDT.  Insufficient  concentration  measurements  were  available  beyond  2215  MDT 
to  evaluate  the  turbulence  burst  effects  on  the  distribution  of  in-plvime 
concentrations . 
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burst  during  Trial  P142154. 


The  velocity  and  concentration  data  taken  during  Trial  P142154  are  incon¬ 
sistent  with  the  intermittent  nocturnal  turbulence  event  model  described  by 
Businger  (1973).  An  Rp-triggered  turbulence  event  in  a  vertically  stratified, 
horizontally  homogeneous  atmosphere  would  have  produced  a  well -mixed  plume 
concentration  field  at  the  tower  beginning  at  2215  MDT.  The  large  horizontal 
velocity  component  perturbations  accompanying  the  disappearance  of  propylene 
at  the  tower  suggests  a  bulk  horizontal  and/or  vertical  displacement  of  the 
plume  rather  than  a  uniform  mixing  event. 

Figure  15  shows  a  second  set  of  time  series  plots  of  wind  components  and 
speed  of  sound  for  two  consecutive  10 -min  time  blocks  which  illustrate  a 
turbulence  burst  observed  during  Trial  P150010.  Table  7  also  gives  the 
statistical  summaries  for  time  blocks  2150013  and  2150023.  The  initial 
turbulence  levels  during  Trial  P150010  were  much  higher  than  those  observed 
during  Trial  P142154,  and  the  stratification  was  not  sufficient  to  reach  the 
critical  Rp.  Nevertheless,  a  turbulence  burst  characterized  by  sudden  wind 
field  perturbations  and  large  changes  in  speed  of  sound  occurred,  as  shown  in 
the  second  10-min  time  block  of  Figure  15.  This  event  triggered  a  DR  increase 
from  0.07,  indicating  virtually  no  contribution  from  convection- dominated 
quadrants ,  to  a  value  of  0.46.  The  turbulence  level ,  as  indicated  by  the  TDI , 
also  increased  by  a  factor  of  four.  The  Businger  (1973)  turbulence  burst 
model  once  again  fails  to  fit  the  observed  scenario , 

An  alternate  explanation  of  the  sporadic  noturnal  turbulence  events 
observed  over  the  high  desert  is  the  collision  between  pools  of  air  within  the 
nocturnal  inversion  layer.  The  high  desert  nocturnal  boundary  layer  contains 
shallow  pools  of  air  with  varying  densities  that  form  over  source  regions  with 
non-uniform  cooling  characteristics.  These  pools  drift  downslope  or  are 
advected  by  the  s3moptic  flow  and  eventually  collide  with  one  another;  the 
distinguishing  characteristics  of  these  collisions  are  the  abrupt  temperature 
and  velocity  component  perturbations  that  occur  during  the  lateral  and 
vertical  displacements  of  the  warmer  pool  by  the  cooler  one.  The  principal 
effect  of  these  collisions  on  nocturnal  dispersion  is  the  bulk  shearing  and/or 
displacement  of  any  puff  or  plume  material  located  at  pool  boundaries.  Also, 
the  turbulence  generated  during  collisons  dissipates  in  an  intermittent 
cascade  that  temporarily  accelerates  in-plume  diffusion. 

5.4  TURBULENT  MIXING  AND  CONCENTRATION  PDFS 

The  effects  of  turbulent  mixing  on  a  diffusing  cloud  are  conveniently 
illustrated  by  PDF  plots.  A  plot  of  an  instantaneous  point-source  concentra¬ 
tion  PDF  displays  the  salient  time -independent  attributes  of  the  concentration 
time  series.  PDF  plots  from  measurements  obtained  under  a  range  of  meteoro¬ 
logical  conditions  at  several  different  positions  within  dispersing  plumes 
illustrate  the  height-  and  stability- dependent  effects  of  turbulent  dispersion 
on  the  distribution  of  plume  concentrations.  Yee  et  al.  (1994)  discuss  the 
various  shapes  that  PDF  curves  can  take,  depending  on  the  exposure  of  a 
dispersing  plume  to  the  different  turbulence  scales  and  intensities.  For 
example,  exponential  decay  from  a  prominent  peak  at  zero  concentration  and  an 
extended  tail  indicate  a  poorly-mixed  plume  exposed  to  weak  small-scale 
turbulence.  Conversely,  strong  mixing  produces  a  non-zero  "shoulder"  or 
secondary  peak  around  x/C  of  one  (where  x  is  the  instantaneous  concentration 
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Horizontal  (u  and  v)  1.5-m  wind  components,  the  3.0-m  vertical  (w)  wind  components,  and  the 
3.0-m  speed  of  sound  (C)  for  two  successive  10-min  periods  illustrating  fully  developed 
nocturnal  turbulence  followed  by  a  turbulence  burst  during  Trial  P150010. 


and  C  is  the  mean  concentration) ,  with  a  truncated  tail  due  to  elimination  of 
high  internal  plume  concentrations.  As  dispersing  material  becomes  thoroughly 
mixed,  the  PDF  assumes  the  shape  of  a  Gaussian  distribution  centered  on  a  x/C 
of  one . 

The  PDFs  presented  below  were  derived  using  a  procedure  described  by  Yee 
et  al.  (1994).  This  procedure  involves  sorting  the  normalized  concentrations, 
X/C  into ^histogram  bins.  Bin  counts  are  then  normalized  by  bin  interval  and 
sample  size.  The  result  is  a  PDF  of  normalized  concentration  frequency  f(x/C) 
plotted  against  x/C,  as  in  Figure  16.  Each  plot  contains  a  set  of  PDFs  from 
TIP-SJ2  detector  measurements  made  at  the  indicated  heights  (0.5  to  3.0  m 
AGL).  Unfortunately,  the  number  of  PDF  plots  available  for  analysis  is  small. 
Four  of  the  available  plots  plus  accompanying  wind  and  turbulence  data  are 
discussed  below  as  case  studies. 

The  PDF  plots  presented  in  Figure  16  are  from  Trials  P070435  (Figure  16a) 
and  P030745  (Figure  16b).  Figure  16a  shows  an  exponential  decrease  in 
concentration  at  all  levels,  with  the  long  "tails"  of  high  concentrations 
characteristic  of  a  relatively  unmixed  plume.  The  tails  are  shortest  at  the 
lowest  level  (0.5  m)  and  increase  monotonically  with  height,  reflecting  the 
effect  of  near -surface  mixing,  even  though  the  measurements  were  made  at  a 
short  (25-m)  travel  distance  from  the  source.  Table  8  provides 
micrometeorological  summaries  from  two  adjacent  10-min  time  blocks  (2070437 
and  2070447)  within  Trial  P070435,  Turbulence  levels  within  these  time 
blocks,  as  indicated  by  the  variances  and  TDI,  are  fairly  typical  for  the  high 
desert  nocturnal  boundary  layer.  The  Obukhov  length  (calculated  using 
rather  than  u.)  and  DR  indicate  stable  nocturnal  conditions  with  little  change 
from  one  time  block  to  the  other. 

In  contrast  to  Figure  16a,  the  PDFs  obtained  from  concentration  detectors 
nearest  the  surface  on  Trial  P030743  (Figure  16b)  show  a  "shoulder"  developing 
in  the  concentration  PDFs  for  the  lower  measurement  levels.  This  shoulder  has 
a  x/C  secondary  maximum  in  the  range  of  1  to  6,  with  rapidly  diminishing 
frequency  at  higher  x/C.  PDFs  from  detectors  positioned  at  higher  levels  have 
longer  tails  and  lack  the  shoulder.  This  distribution  of  PDFs  occurred  during 
diffusion  over  a  distance  of  50  m  in  a  boundary  layer  in  transition,  with 
surface  heating  inducing  incipient  convective  turbulence  at  the  lowest  levels 
and  the  higher  levels  remaining  stable.  The  corresponding  micrometeorological 
summaries  in  Table  8  include  a  preconvection  10-min  block  (2030743)  followed 
by  a  10-min  block  during  which  convection  penetrated  to  the  1.5-m  level 
(2030753).  The  change  in  1.5-m  DR  from  0.35  (stable,  but  approaching  neutral) 
in  the  first  10-min  block  to  2.29  (convection- dominated)  in  the  next  10-min 
block  clearly  defines  a  rapid  stability  progression  during  the  transition. 

The  stability  transition  at  the  6-m  level  is  delayed  in  time,  with  a  DR  of 
0.02  (strongly  stable)  in  the  first  time  block  increasing  to  0.41  (approaching 
neutral)  in  the  second.  Although  PDFs  from  detectors  mounted  at  6  m  are  not 
available,  it  is  evident  from  the  2.5-  and  3.0-m  PDF  plots  in  Figure  16b  that 
surface -based  convection  was  not  penetrating  above  the  first  2  m.  The  3-m 
Obukhov  length  L  is  negative  for  both  time  blocks,  indicating  that  the  L- 
derived  transition  from  stable  to  unstable  occurred  somewhat  earlier  than  the 
DR- derived  transition. 
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Figure  17  shows  a  second  set  of  concentration  PDFs.  These  PDFs  are  for 
concentration  profiles  taken  during  Trials  P142154  (Figure  17a)  and  P150010 

17b) .  Trial  P142154  is  representative  of  diffusion  over  a  distance  of 
100  m  into  a  quiescent  nocturnal  flow  that  was  interrupted  by  a  turbulence 
burst  in  mid- trial  (see  Figure  14).  The  PDFs  in  Figure  17a,  with  their 
exponential  slope  and  long  tail  to  high  x/C,  are  similar  in  appearance  to  the 
nocturnal  case  presented  in  Figure  16a  except  for  the  flatness  of  the  PDF  at 
the  4.0-m  level  in  Figure  17a.  The  flatness  of  the  4.0-m  PDF  for  x/C  of  7  to 
30  Indicates  very  little  internal  mixing,  which  is  characteristic  of  a  plume 
within  a  quiescent  layer .  Plume  material  was  present  on  the  tower 
mainly  during  the  first  part  of  the  trial  (represented  by  the  2142205  time 
period  block  in  Table  7),  prior  to  the  turbulence  burst.  The  air  pool  colli¬ 
sion  that  produced  the  turbulence  also  displaced  the  plume,  with  only  very 
intermittent  weak  concentrations  occurring  at  the  tower  for  the  remainder  of 
the  trial.  Consequently,  the  PDFs  in  Figure  17a  are  representative  only  of 
diffusion  within  the  quiescent,  pre-burst  turbulence  field. 

The  micrometeorological  svimmaries  shown  in  Table  7  for  the  2142205  and 
2142215  time  blocks  document  changes  that  occurred  within  a  quiescent  noctur¬ 
nal  ABL  during  a  turbulence  burst.  The  velocity  and  temperature  variances 
increase  dramatically,  reflected  by  a  proportionate  TDI  increase.  The  1.5-m 
level  DR  indicates  a  slight  trend  from  stable  (0.10)  toward  neutral  (0.33) 
during  this  time.  In  contrast,  the  change  in  Obukhov  length  at  3.0-m  from 
-0.6  to  +1.8  does  not  accurately  reflect  stability  conditions  or  the  occur¬ 
rence  of  a  stable  case  turbulence  burst.  Again,  it  is  not  possible  to  compare 
from  pre-  and  post- turbulence  burst  conditions  because  the  plume  was 
displaced  away  from  the  tower  as  a  consequence  of  this  event. 

In  contrast  to  Figure  17a,  Figure  17b  (Trial  P150010)  shows  PDFs  for 
diffusion  over  a  100  m  distance  into  a  stable  boundary  layer  with  initial 
turbulence  levels  similar  to  those  of  Trial  P070435.  Turbulence  levels  for 
Trial  P150010  are  represented  by  the  initial  10-min  period  beginning  at 
2150013,  followed  by  an  abrupt  turbulence  increase  at  2150023  (see  Table  7  and 
15) .  As  in  Trial  P142154,  the  turbulence  event  displaced  the  plume 
away  from  the  tower .  Consequently ,  the  available  PDFs  are  most  representative 
of  diffusion  into  the  pre-event  nocturnal  ABL.  Turbulence  levels  during  this 
trial  were  relatively  high  from  the  beginning,  as  indicated  by  the  well  mixed 
PDFs  at  the  0.5- ,  1.0- ,  and  2.0-m  levels.  Surface -generated  turbulence 
apparently  diminished  rapidly  with  height,  leaving  a  4-m  PDF  indicative  of  a 
relatively  poorly  mixed  plume. 

The  Trial  P150010  turbulence  event  beginning  at  2150023  produced  turbu¬ 
lence  statistics  comparable  to  the  statistics  obtained  for  the  turbulence 
event  of  Trial  P142154.  The  Trial  P150010  1.5-m  DR  responded  in  a  similar 
fashion  to  that  of  Trial  P142154  by  increasing  from  0.02  to  0.44  as  a  conse¬ 
quence  of  the  turbulence  event.  The  Obukhov  length  again  provided  no  useful 
indication  of  the  increase  in  turbulence  between  the  2150013  and  2150023  time 
periods . 
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SECTION  6.  CONCLUSIONS  AND  RECOMMENDATIONS 


Improved  measurement  and  analysis  techniques  are  beginning  to  reveal  the 
detailed  physics  of  atmospheric  dispersion  processes  that  heretofore  have  been 
subject  to  guesswork  and  gross  parameterizations .  In  particular,  the  study  of 
dispersion  near  surfaces  has  been  hindered  by  the  absence  of  measurement 
techniques  able  to  resolve  small-scale  processes  in  an  environment  character¬ 
ized  by  strong  gradients  and  non- steady  fluxes  of  heat  and  momentum.  This 
region  is  currently  beyond  the  reach  of  numerical  modeling  techniques  that  are 
successfully  resolving  larger  scale  atmospheric  motions.  The  Evaporation, 
Recirculation,  and  Dispersion  in  Light  Winds  Dispersion  test:  (1)  generated 
data  bases  for  validating  evaporation  and  recirculation  models,  and  (2) 
provided  a  first  look  into  the  physics  of  some  near -surface  high  Reynolds 
number  processes .  Much  work  remains  to  achieve  a  clear  understanding  of  these 
processes . 

The  Evaporation  Subtest  provided  a  valuable  data  set  detailing  the  effects 
of  droplet  size,  thickener  concentration,  and  surface  porosity  on  the  evapora¬ 
tion  of  the  chemical  agent  simulants  DEM  and  MES.  A  detailed  description  of 
the  Evaporation  Subtest  and  an  evaluation  of  subsequent  modeling  results  is 
provided  by  Roberts  (1996).  Roberts'  evaluation  of  an  analytical  droplet 
evaporation  prediction  model  found  good  initial  results ,  but  a  general 
evaporation  rate  overprediction  with  increasing  evaporation  time. 

The  Recirculation  Subtest  characterized  flow  around  a  simple  bluff  body 
and  the  dissipation  of  windborne  tracer  gas  in  its  wake.  Existing  recircula¬ 
tion  region  contaminant  accumulation  models  were  evaluated.  Principal 
conclusions  are:  (1)  wake  width  is  a  more  pertinent  dimension  than  obstacle 
width,  and  (2)  stability  has  minimal  effect  on  the  properly  non- dimens  ion - 
alized  residence  time.  Mavoridis  (1997)  uses  these  and  related  study  results 
to  develop  improved  urban  dispersion  model  algorithms. 

Analyses  of  turbulence  and  dispersing  gas  concentration  profiles  within 
the  buffer  and  inertial  sublayers  have  revealed  some  of  the  physics  of 
important  non- steady  processes  such  as  the  transition  through  sunrise  and  the 
nocturnal  turbulence  burst.  Fourier  and  u'r'  covariance  quadrant  analyses  were 
found  to  provide  useful  insights  into  these  non-steady  processes.  The 
transition  through  sunrise  was  examined  in  detail  to  determine  how  the 
atmosphere  evolves  from  a  stable  nocturnal  state  to  a  state  of  convection. 

The  adiabatic  condition,  with  no  vertical  exchange  of  heat,  often  has  been 
presumed  to  exist  during  this  transition.  However,  rather  than  existing  in  a 
true  adiabatic  state,  the  desert  atmosphere  in  transition  occurs  when  there  is 
a  net  balance  between  eddies  transporting  heat  to  and  from  the  surface.  An 
intense  near- surface  convective  layer  develops  after  sunrise  and  works  against 
the  trapping  inversion  lid.  This  depiction  of  boundary  layer  transition 
differs  profoundly  from  the  well  mixed  adiabatic  model  that  pertains  to 
regions  like  Northern  Europe.  The  desert  transition  pattern  implies  that 
vigorous  mixing  can  occur  within  confined  layers,  creating  transient  episodes 
of  high  ground- level  pollutant  concentrations  from  plumes  trapped  in  elevated 
stable  layers.  Similar  analyses  of  the  sunset  transition  are  needed  to 
complete  the  description  of  these  non- steady  processes. 
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The  lack  of  satisfactory  similarity  characterization  of  the  non- steady 
lower  portion  of  the  surface  boundary  layer  led  to  the  development  of  two  new 
indices,  the  Diabatic  Ratio  (DR)  and  the  Total  Diabatic  Influence  (TDI) .  DR, 
the  ratio  of  convective  to  shear-driven  turbulence,  is  a  stability  indicator, 
while  TDI  provides  an  estimate  of  the  magnitude  of  vertical  turbulent  motions. 
Both  indices  are  needed  to  adequately  describe  the  turbulent  state  of  the 
surface  boundary  layer. 

DR  and  TDI  appear  to  be  more  robust  indicators  than  Obukhov  length  L, 
particularly  with  light  winds  and  non- steady  conditions  within  the  lower 
portion  of  the  surface  boundary  layer.  An  advantage  offered  by  the  methods 
used  to  compute  DR  and  TDI  is  that  they  do  not  involve  the  summing  of  positive 
and  negative  quantities  required  for  flux  computations.  The  result  of 
algebraically  summing  positive  and  negative  covariance  contributions  is  a  flxix 
that  is  often  small  in  comparison  with  its  uncertainty.  Using  a  ratio  (the 
DR)  or  an  absolute  sum  (the  TDI)  of  covariance  quadrant  contributions  creates 
indices  that  are  relatively  large  and  tractable,  while  remaining  physically 
relevant.  These  indices  offer  the  added  advantage  that  they  can  be  easily 
computed  in  an  operational  setting  using  near-surface  measurements  of  the 
horizontal  wind  and  temperature  components  from  a  2 -axis  sonic  anemometer/ 
thermometer.  DR  and  TDI  appear  to  be  statistically  stable  when  computed  over 
averaging  periods  on  the  order  of  10  min.  Thus,  they  could  provide  near  real¬ 
time  estimates  of  the  turbulent  state  of  the  surface  boundary  layer  in  an 
operational  setting.  These  proposed  indices  are  not  presently  supported  by 
boundary  layer  theory,  and  they  have  not  been  subject  to  independent  testing. 
Further  investigations ,  to  include  development  of  a  supporting  theory  and 
applications  testing,  are  needed  to  determine  their  overall  usefulness  as 
surface  layer  stability  and  turbulence  indicators. 
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SECTION  7 ■  APPENDICES 


APPENDIX  A.  CONSOLIDATED  TIP-SJ  AND 
MICROMETEOROLOGICAL  DATA 


The  statistical  stunmaries  presented  in  Appendix  A  were  obtained  from  sonic 
anemometer/thermometers  and  TIPSJ2  photoionization  detectors  mounted  on  a  9-m 
tower  at  the  indicated  distances  downwind  from  a  continuous  point  source  of 
propylene.  The  trial  name  convention  is  JJJHHMM,  where  JJJ  is  the  Julian  date 
and  HHMM  is  the  nominal  trial  start  time  (MDT) .  The  actual  trial  start  and 
end  times  are  given  in  Universal  Coordinated  Time  (UTC) .  Wind  component, 
sonic  temperature  statistics,  and  temperature  profiles  are  presented  in  the 
top  half  of  each  summary  page.  The  covariance  quadrant  table  includes  u'T' 
covariance  components  for  the  1.5-  and  6-m  levels,  and  (7F  covariance  compo¬ 
nents  for  the  3-  and  9-m  levels.  With  the  exception  of  skewness  and  kurtosis, 
all  of  the  u  and  v  statistics  were  computed  after  coordinate  rotation  into  the 
mean  wind.  Below  the  micrometeorological  summary  is  a  profile  plot  of  peak 
and  mean  concentrations,  normalized  by  the  highest  profile  mean  concentration. 
Concentration  profile  data  are  tabulated  to  the  right  of  the  plot.  Trial 
concentration  data  were  provided  by  Dr.  Eugene  Yee  of  DRES. 
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Trial:  P030545  Micrometeorological  and  Concentration  Data  Summary  Date:  22  July  1995 

Source  Distance:  50  m  Source  Height:  1.0  m  Peak  Concentration  Ratio; 

Time  (UTC) :  1146  to  1159  Obukhov  Length  (m) :  +2.  Temperature  (°C)  0-m:  12.9  0.5-m:  13.3  2.0-m:  13. 
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Trial;  P030745  Micrometeorologlcal  and  Concentration  Data  Summary  Date:  22  July  1995 

Source  Distance:  50  m  Source  Height:  0.05  m  Peak  Concentration  Ratio:  7.94 
Time  (UTC) :  1343  to  1413  Obukhov  Length  (m) :  -9.  Temperature  (“C)  0-m:  17.9  0.5-m:  20.2  2.0-m:  19.2 
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Trial.  F060530  Micrometeorological  and  Concentration  Data  Summary  Date:  25  July  1995 

Source  Distance:  25  m  Source  Height:  0.05  m  Peak  Concentration  Ratio:  ND 
Time  (UTC) :  1128  to  1143  Obukhov  Length  (m) :  +0.1  Temperature  (°C)  0-m:  10.8  0.5-m:  10.9  2.0-m:  11.8 
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Concentrations 


Trial:  P070435  Micrometeorological  and  Concentration  Data  Sunraiary  Date:  26  July  1995 

Source  Distance:  25  m  Source  Height:  1.0  m  Peak  Concentration  Ratio:  9.94 
Time  (UTC) :  1032  to  1107  Obukhov  Length  (m) :  12.  Temperature  (°C)  0-m:  16.4  0.5-m:  19.8  2.0-m:  20. 
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Trial:  P070605  Micrometeorological  and  Concentration  Data  Summary  Date:  26  July  1995 

Source  Distance:  20  m  Source  Height:  0.05  m  Peak  Concentration  Ratio:  ND* 

Time  (UTC) :  1200  to  1235  Obukhov  Length  (m) :  +10.  Temperature  (°C)  0-m:  13.9  0.5-m:  15.7  2.0-m: 


o 

o 

in 

o 

in 

O 

in 

in 

CM 

m 

CM 

CM 

rH 

fH 

o 

o 

A-13 


(Insufficient  Data) 


A- 14 


00 


I 


0)  (O  in  ^ 

^  o  OJ 

d  d  d  d 


!  r  i  ® 
^  d 


^  ^ 
in  5  to  Q 

6^0^ 


S  Q  S  D 

2  2  d  2 


O) 

CO 

CO 

o 

CO 

CM 

CO 

CO 

CO 

CM 

d 

d 

d 

d 

CD 

CO 

1 

O 

, 

•— 

1 

O 

• 

o 

d 

d 

in 

in 

O 

o 

r- 

o 

w“ 

z 

z 

d 

o 

o 

in 

o 

CD 

D 

z 

Z 

o 

d 

a  “j  Q 

in  Z  o  2 


§  g  S  g 

CN  ^  CM  ^ 


^  Q  ®  Q 

t-'  Z  Z 


in  q  q  q 
^  CO  CO*  cn 


A- 16 


Trial:  P100630  Micrometeorological  and  Concentration  Data  Suiranary  Date:  29  July  1995 

Source  Distance:  75  m  Source  Height:  0.05  m  Peak  Concentration  Ratio:  ND* 

Time  (UTC)  :  1226  to  1302  Obukhov  Length  (m) :  Temperature  (°C)  0-m:  14.5  0.5-m:  15.6  2.0-m:  16.5 
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APPENDIX  C.  LIST  OF  SYMBOLS 


c  Speed  of  sound 
C  Mean  measured  concentration 

Cuw  Correlation  between  horizontal  and  vertical  velocity  fluctuations 
Correlation  between  temperature  and  vertical  velocity  fluctuations 
Cut  Correlation  between  temperature  and  horizontal  velocity  fluctuations 
Cp  Specific  heat  at  constant  pressure 
CpK  Peak  (99th  percentile)  measured  concentration 
d  Droplet  diameter 

D21  Horizontal  component  of  turbulent  thermal  diffusivity 
D22  Vertical  component  of  turbulent  thermal  diffusivity 
g  Gravitational  acceleration 

H  Heat  Flux, 

h  Cube  face  length 
h'  Wake  width 

h^  Momentum  deficit  thickness 
hr  Surface  texture  or  grain  size 

i  Concentration  fluctuation  intensity  (ratio  of  the  concentration 
standard  deviation  to  the  mean  concentration) 

K  Kurtosis  of  the  concentration  field  measured  at  a  detector 

k  von  Karman  constant  (nondimens ional  parameter  with  a  magnitude  of  0.4) 

L  Obukhov  length  (L  pCpU*^/kgH) 

n  Sample  size 

P/M  Ratio  of  plume  profile  peak  concentration  to  plume  centroid  mean  concen¬ 
tration 

Turbulent  kinetic  energy 
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LIST  OF  SYMBOLS 
(CONTINUED^ 

Rp  Flux  Richardson  number 
Rp  Plume  Richardson  Number 

R^  Momentum  deficit  thickness  Reynolds  number 

S  Skewness  of  the  concentration  distribution  measured  at  a  detector 

T  Temperature,  degrees  Kelvin 

T '  ^  Temperature  variance 

t  Averaging  time 

t^  Residence  time 

t.  Nondimens ional  residence  time 

u  Wind  speed 

U  Mean  wind  speed  at  plume  centroid 
Uy  Surface  shear  stress 

u*  Friction  velocity,  the  square  root  of  the  vertical  momentum  flux 
(-uV) 

u'  Axial  velocity  root-mean  square  (rms) 
u'^  Alongwind  velocity  variance 

U„  Free  stream  velocity  above  the  boundary  layer 
u*  Inner-normalized  mean  axial  velocity  (u/u,.) 
u*  Inner-normalized  axial  velocity  rms  (u'/uj 
v'^  Crosswind  velocity  variance 
w  Vertical  velocity 
w'^  Vertical  velocity  variance 
Wp  Free  convection  scale 
w.  Convective  velocity  scale 

X  Downwind  distance  from  dissemination  point  to  detector  array 
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LIST  OF  SYMBOLS 
(CONTINUED^ 


X  Dimensionless  downwind  distance  in  a  convective  boundary  layer 


y"^  Spanwise  pocket  scale  (yuy^) 

Zc  Height  of  the  plume  centroid 

Zi  Mixed  Layer  depth,  height  of  the  atmospheric  boundary  layer 

Zo  Aerod3niamic  roughness  parameter 

Zp  Height  of  the  plume  peak  concentration 

Zr  Height  of  detector  above  ground 

Zs  Height  of  source  above  ground 

z*  Nondimens ional  height 

€  Turbulent  kinetic  energy  dissipation  rate 
Xj  Inertial  scale 
Xt  Taylor  microscale 

6  Potential  temperature,  degrees  Kelvin 

T  Plume  travel  time  normalized  by  the  Lagrangian  timescale 
p  Density  of  air 

V  Tracer  particle  vertical  drift  velocity 

V  Viscosity  of  air 
Vj  Eddy  viscosity 

X  Instantaneous  concentration 
i?  Mean  tracer  particle  drift  velocity 
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HEADQUARTERS 

AFSOC/CEX 

100  BARTLEY  STREET,  SUTIE  218-E 
ATTN;  MR.  TOM  GRAHAM 
HULBERT  AFB,  FL  32544-5273 

HQ  USCENTCOM 
CCJ3-OG 

ATTN:  MAJ  RAY  VAN  PELT 
MACDILL  AFB,  FL  33594 

•  •  •  * 


COMMANDER  NAVAL  AIR  FORCE 
U.S.  PACinC  FLEET 
BOX  357051 

SAN  DIEGO,  CA  92135-7051 

COMMANDER,  NAVAL  SURFACE  FORCE 
U,S.  PACinC  FLEET  (ATTN;  N66) 

NAVAL  AMPHIBIOUS  BASE,  CORONADO 
SAN  DIEGO,  CA  92155-5035 

COMMANDER  SEVENTH  FLEET 
FPO  AP  96601 

COMMANDER  THIRD  FLEET 
FPO  AP  96601-6001 

COMMANDER  TRAINING  COMMAND 
U.S.  PACinC  FLEET 
SAN  DIEGO,  CA  92147 

COMMANDER  TRAINING  COMMAND 
U.S.  ATLANTIC  FLEET 
ATTN;  N32 
NORFOLK,  VA  23511 

COMMANDER 

U.S.  ATLANTIC  FLEET  (N7) 

1562  MITSCHER  AVE,  SUITE  250 
NORFOLK,  VA  23551-2487 

DEPARTMENT  OF  THE  NAVY 
FLEET  TRAINING  GROUP 
DETACHMENT  NORFOLK 
1310  8TH  ST. 

NORFOLK,  VA.  23521-5150 

COMMANDER  NAVAL  SURFACE  FORCE 
US  ATLANTIC  FLEET  N541A 
1430  MITSCHER  AVE 
NORFOLK,  VA  23551-2494 


COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
ATTN:  (MR  S.  ENATSKY) 

SEA 0566  (MRP, JUNG) 

2531  JEFFERSON  DAVIS  HWY 
WASHINGTON,  DC  22242-5160 

COMMANDER 

NAVAL  SPECIAL  WARFARE  COMMAND 

2000  TRIDENT  WAY 

SAN  DIEGO,  CA  92155-5000 

NAVAL  SURFACE  WARFARE  CENTER 
CODE  B-50  (MR.  J.  HORTON) 

CODE  B-53  OHR.  BYRNE) 

CODE  B-51  (MR.  YENCHA) 

CODE  B-S3  (MR.  POMPEII) 

DAHLGREN,  VA  22448-5000 

COMMANDING  OFFICER 

NAVAL  TRAINING  SYSTEMS  CENTER 

CODE  4972  (T.  FRANZ) 

12350  RESEARCH  PARKWAY 
ORLANDO,  FL  32826-3224 

DIRECTOR,  ONI 
ATTN:  ONI  2323 
4251  SUITLAND  ROAD 
WASHINGTON  D.C.  20395-5720 

COMMANDER 

NAVAL  HEALTH  RESEARCH  CENTER 
ATTN:  DR  ROBERT  PAZOS 
BOX  85122 

SAN  DIEGO  CA  92138 

OFFICER  IN  CHARGE 

NAVAL  RESEARCH  LAB  DETACHMENT 

CBD  COM 

ABERDEEN  PROVING  GROUND,  MD  21010-5323 

BUREAU  OF  MEDICINE  AND  SURGERY 
ATTN:  MED-OZC 
2300  E  STREET,  NW 
WASHINGTON,  DC  20372-5300 

COMMANDING  OFFICER 
NAVAL  RESEARCH  LABORATORY 
CODE  6116  R  A  LAMONTAGNE 
WASHINGTON,  DC  20375-5320 

NAVAL  AIR  WARFARE  CENTER-AIRCRAFT  DIVISION 

ATTN:  GAITHER  M  PUGH  4.1.8 

UNIT  5  BLDG  2187  SUITE  1180 

48110  SHAW  ROAD 

PATUXENT  RIVER  MD  20670-1906 

SURFACE  WARFARE  DIVISION 
ATTN:  CODE  N862,  RM  4D537 
WASHINGTON,  D.C.  20350-2000 
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UTC  (SCW)  ECKHOFF 
5510  CBC  8TH  STREET 
GULFPORT,  MS  39501-5003 


•  •  •  • 

COMMANDANT 
U.S.  MARINE  CORPS 
ATTN:  APW  R-P 
WASHINGTON,  DC  20380 

COMMANDING  OFRCER 

MARINE  AIRCRAFT  GROUP  31  (S-3/NBC) 

2D  MARINE  AIRCRAFT  WING,  FMF,  ATLANTIC 
BEAUFORT,  SC  29904-0120 

COMMANDING  OFFICER 

MARINE  AIR  CONTROL  GROUP  28  (S-3/NBC) 

2D  MARINE  AIRCRAFT  WING,  FMF  ATLANTIC 

P.O.BOX  8069 

CHERRY  POINT,  NC  28533 

COMMANDING  OFFICER 
MARINE  WING  SUPPORT  GROUP  27  (S-3/NBC 
2D  MARINE  AIRCRAFT  WING,  FMF  ATLANTIC 
CHERRY  POINT,  NC  28533 

COMMANDING  OFFICER 
MARINE  AIR  WING  HEADQUARTERS 
SQUADRON  2  (S-3/NBC) 

2D  MARINE  AIRCRAFT  WING,  FMF,  ATLANTIC 
CHERRY  POINT,  NC  28533 

COMMANDING  GENERAL 

MARINE  CORPS  COMBAT  DEVELOPMENT  COMMAND 
ATTN:  (CWO  TURVILLE)  C  443 
2042  BROADWAY  ST,  SUITE  17 
QUANTICO,  VA  22134-5021 

COMMANDING  GENERAL 
FLEET  MARINE  FORCES  PACIFIC 
ATTN:  G3/NBC 
CAMP  HJVl.  SMITH,  HI  96861 

COMMANDING  GENERAL 
FLEET  MARINE  FORCES  ATLANTIC 
ATTN:  G3/NBC 
1468  INGRAM  STREET 
NORFOLK,  VA  23551-2596 

USMC  LIAISON  OFFICER  TO  TECOM 
ATTN:  AMSTE-ML 

ABERDEEN  PROVING  GROUND,  MD  21005 

COMMANDER 
MARCORSYSCOM 
ATTN:  AW 

2033  BARNETT  AVENUE,  SUITE  315 
QUANTICO,  VA  22134-5010 

COMMANDING  GENERAL 
II  MARINE  EXPEDITIONARY  FORCE 
ATTN:  G-3/NBC,  PSC  BOX  20080 
CAMPLEJEUNE,NC  28542-5401 


COMMANDING  OFFICER 

MARINE  WING  SUPPORT  SQUADRON  271  (S-S/NBC) 
2D  MARINE  AIRCRAFT  WING,  FMF,  ATLANTIC 
MCAS  CHERRY  POINT,  NC  28533 

COMMANDING  OFFICER 

MARINE  WING  SUPPORT  SQUADRON  274  (S-3/NBC) 
2D  MARINE  AIRCRAFT  WING,  FMF,  ATLANTIC 
MCAS  CHERRY  POINT,  NC  28533 

COMMANDING  OFFICER 

MARINE  WING  SUPPORT  SQUADRON  273  (S-3/NBC) 
2D  MARINE  AIRCRAFT  WING,  FMF,  ATLANTIC 
MCAS  BEAUFORT,  SC  29904-6035 

COMMANDING  OFFICER 
HQ,  MAG-26,  (NBC) 

PSC  BOX  21010 

JACKSONVILLE,  NC  28545-1010 

COMMANDING  GENERAL 
ATTN  NBC  DEFENSE  OFFICER 
III  MARINE  EXPEDITIONARY  ORCE 
UNIT  35601 
FPO  AP  96606-5601 

COMMANDING  OFFICER 

31ST  MARINE  EXPEDITIONARY  UNIT  FMF 

ATTN  NBCO  CW03  TERRY  MURPHY 

UNIT  35621 

FPOAP  96606-5621 

UNITED  STATES  MARINE  CORPS 
2ND  MARINE  AIRCRAFT  WING  (G3/NBC) 

U.S.  MARINE  CORPS  FORCES,  ATLANTIC 
POSTAL  SERVICES  BOX  8050 
CHERRY  POINT,  NC  28533-0050 

UNITED  STATES  MARINE  CORPS 
MARPORLANT  NBC 
CAMPLEJEUNE,  NC  28542-5701 

NBC  OFFICER 
MWHS-1 1ST  MAW 
UNIT  37121  (G-3/NBC) 

FPO  AP  96603-7121 

•  •  *  • 


ADMINISTRATOR,  DEFENSE  TECHNICAL 
INFORMATION  CENTER  (2  copies) 

CAMERON  STATION 

8725  JOHN  J.  KINGMAN  RD  SUITE  0944 

FORT  BEL  VOIR,  VA  22060-6218 

DIRECTOR 

CENTERS  FOR  DISEASE  CONTROL 
ATTN:  LOGGING  CONTROL  OFFICER  -  1/B413 
(for  ASST  DIR,  EDWIN  K.  GRAY) 

ATLANTA,  GA  30333 


COMMANDING  GENERAL 
2D  MARINE  DIV  FMF 
PSC  BOX  20003 
CAMPLEJEUNE,NC  28540 
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